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LEVEE  UNDERSEEPAGE  ANALYSIS  FOR  SPECIAL 
FOUNDATION  CONDITIONS 


PART  I:  INTRODUCTION 

Background,  Purpose,  and  Scope 

1.  A  Repair,  Evaluation,  Maintenance,  and  Rehabilitation  (REMR)  Levee 
Underseepage  Workshop  was  held  at  the  US  Army  Engineer  Waterways  Experiment 
Station  (WES)  on  19  April  1984  to  establish  research  needs  related  to  control 
of  levee  underseepage.  Representatives  from  the  Rock  Island,  St.  Louis, 
Memphis,  and  Vicksburg  Corps  of  Engineers  (CE)  Districts  attended  the  work¬ 
shop.  One  of  the  research  tasks  established  was  comparing  predicted  levee 
underseepage  conditions  to  observed  performance.  In  September  1986,  a  criti¬ 
cal  review  of  ur.Jerseepage  analysis  procedures  was  prepared  by  this  author 
(Wolff  1986)  under  an  Interagency  Personnel  Agreement  with  the  WES  Geotechni¬ 
cal  Laboratory.  The  workshop  and  the  critical  review  both  indicated  that 
levee  and  foundation  geometry  may  significantly  affect  seepage  conditions,  and 
the  method  of  modeling  geometry  in  analysis  may  affect  performance 
predictions . 

2.  The  focus  of  underseepage  analysis  is  to  estimate  or  predict  the 
residual  head,  h^  ,  and  exit  gradient,  i  ,  at  the  landside  levee  toe  or  berm 
toe  during  high  water.  Where  the  predicted  gradient  exceeds  some  critical 
value,  typically  0.85,  control  measures  such  as  relief  wells  or  seepage  berms 
are  provided.  The  analysis  procedure  presently  used  by  the  CE  WES  1956a; 

US  Army,  Office,  Chief  of  Engineers  1978)  assumes  two-dimensional  (2-D)  flow 
and  models  the  subsurface  profile  as  two  horizontal  layers  of  uniform  thick¬ 
ness.  In  the  Mississippi  River  Valley,  the  upper  layer,  or  semipervious  top 
stratum,  is  typically  2  to  20  ft*  thick  and  consists  of  clays,  silts,  silty 
sands,  or  combinations  thereof.  The  lower  layer,  or  pervious  substratum,  is 
tvpically  50  to  150  ft  thick  and  consists  of  clean  sand,  gravelly  sand,  or 
siltv  sand.  The  analysis  procedure  is  based  on  a  solution  obtained  by  Bennett 


*  A  table  of  factors  for  converting  nor.-Sl  units  of  measurement  to  SI 
(metric!  units  is  presented  on  page  7. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

degrees  (angle) 
feet 

feet  per  minute 
inches 

miles  (US  statute) 


By 

0.01745329 

0.3048 

0.5080 

2.54 

1.609347 


_ To  Obtain _ 

radians 

metres 

centimetres  per  second 

centimetres 

kilometres 


6.  Use  of  the  developed  computer  programs  is  described  in  Appendices  A, 
B,  and  C;  the  programs  LEVEE3L,  LEVEEIRR,  LEVEECOR  and  the  conventional  method 
are  compared  in  Appendix  D;  the  notation  is  documented  in  Appendix  E. 


Previous  Studies 


Infinitely  long  foundations 

7.  Bennett  (1946)  derived  a  solution  for  steady-state  flow  through  a 
two-layer  foundation  consisting  of  a  semipervious  top  blanket  of  thickness  z 
overlying  a  pervious  substratum  of  thickness  d  ,  both  of  infinite  lateral 
extent.  In  Bennett's  analysis,  flow  is  assumed  downward  vertical  in  the 
riverside  top  blanket,  landward  horizontal  in  the  substratum,  and  upward  ver¬ 
tical  in  the  landside  top  blanket.  Bennett  stated  that  the  substratum  must  be 
at  least  10  times  as  pervious  as  the  top  blanket  for  these  assumptions  to  be 
reasonable;  this  is  almost  always  the  case  for  levees  in  the  Mississippi 
Valley.  Bennett's  analysis  is  summarized  in  Figure  1.  To  calculate  the 

residual  head,  h  ,  at  the  landside  levee  toe  where  the  foundations  layers 
o 

are  of  infinite  length,  the  pervious  substratum  and  semi  pervious  top  blanket 
are  replaced  by  finite  "effective"  lengths  of  pervious  substratum  and  imper¬ 
vious  top  stratum.  These  effective  lengths  are  designated  x^  on  the  river¬ 
side  and  x^  on  the  landside  and  are  a  function  of  the  thicknesses  and 
permeabilities  of  the  two  materials.  The  base  width  of  the  levee  is  desig¬ 
nated  x^  .  With  the  top  stratum  now  impervious,  the  head  in  the  substratum 

varies  linearly  with  distance,  and  the  head  at  the  levee  toe,  h  ,  can  be 

o 

found  by  simple  interpolation: 


h  =  - - -  (1 

O  Xj  +  X^  +  X^ 

where  H  is  the  net  head  on  the  levee  (difference  between  river  stage  and 
landside  ground  or  tailwater) . 

Foundations  of  finite  length 

8.  Solutions  for  effective  entrance  and  exit  distances  for  finite 
foundation  lengths  were  presented  by  Bennett  and  extended  and  summarized  in 
Technical  Memorandum  (TM  3-424  (WES  1964)).  The  semipervious  top  blanket  and 
pervious  substratum  extend  for  distances  of  riverward  of  the  riverside 
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(1946),  which  was  extended  and  summarized  in  WES  Technical  Memorandum  3-424 
(WES  1956a).  It  is  referred  to  herein  as  the  "conventional  analysis." 

Although  actual  foundation  conditions  may  be  highly  irregular,  the  conven¬ 
tional  analysis  requires  the  foundation  geometry  to  be  transformed  into  an 
"equivalent"  system  of  two  horizontal  layers  of  uniform  thickness.  Develop¬ 
ment  of  the  equivalent  foundation  may  involve  changing  both  the  thickness  and 
permeability  of  the  top  blanket.  As  this  transformation  process  can  be  highly 
subjective,  performance  predictions  for  irregular  foundation  conditions  are 
often  difficult  to  make  and  are  unreliable. 

3.  The  "conventional  analysis"  assumptions  allow  a  closed-form  solution 
to  the  differential  equation  for  the  piezometric  head  at  the  interface  between 
the  top  blanket  and  the  substratum.  The  residual  head  and  gradient  could 
readily  be  calculated  using  1950's  techniques  such  as  slide  rules  and  charts. 
Digital  computers  and  numerical  methods  now  allow  solution  of  the  flow 
equation  for  very  complex  conditions.  Finite  element  programs  are  available 
that  can  model  any  conceivable  seepage  problem  (e.g.  Tracy  1973);  however,  they 
are  seldom  used  for  levee  underseepage  analysis  because  the  effort  required 
for  problem  description  and  coding  is  usually  undesirable  for  routine 
calculations.  For  levee  underseepage  problems,  an  analysis  technique  is 
desired  wherein  the  heads  and  gradients  can  be  obtained  as  a  function  of 
relatively  few  parameters,  facilitating  repetitive  calculations  for  numerous 
foundation  sections  along  the  length  of  a  levee. 

4.  The  purpose  of  this  research  was  to  develop  analysis  procedures  that 
are  not  constrained  by  some  of  the  assumptions  in  the  conventional  procedure. 

A  second  purpose  was  to  investigate  whether  improved  performance  predictions 
could  be  made  using  the  developed  procedures.  As  part  of  the  research,  three 
computer  programs  were  developed  to  perform  underseepage  analysis  for  three 
special  but  relatively  common  geometric  conditions.  These  are: 

a.  Program  LEVEE3L  for  analysis  of  foundations  consisting  of  three 
layers  of  uniform  thickness. 

b.  Program  LEVEEIRR  for  analysis  of  foundations  consisting  of  two 
layers  of  irregular  shape  (nonuniform  thickness). 

£.  Program  LEVEECOR  for  analysis  of  underseepage  at  angles  or 
"corners"  in  levee  alignment. 

5.  For  each  of  the  three  types  of  special  foundation  conditions,  two  to 
four  prototype  reaches  were  analyzed  and  the  results  compared  to  actual 
performance  data  (piezometer  readings  during  flood). 
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Dashed  line  ie  piezometric  head  at  base  of  top  blanket 


Figure  i.  Bennett's  analysis 


levee  toe  and  landward  of  the  landside  levee  toe,  respectively,  the 

effective  lengths  x^  and  x^  are  functions  of  and  .  The  pervious 

substratum  may  be  modeled  as  either  open  or  blocked  at  the  entrance  and  exit 
points.  The  analyses  for  open  entrance  and  exit  conditions  are  illustrated  in 
Figure  2. 


Back-Calculation  of  Parameters 


9.  A  number  of  studies  have  been  made  where  residual  heads  and  gradi¬ 
ents  observed  during  high  water  have  been  compared  to  design  calculations 
(e.g.  WES  1964;  US  Army  Engineer  District  (USAED) ,  St.  Louis  1976;  Cunny  1980; 
McClelland  Engineers  1985).  Some  of  the  findings  from  these  studies  have  been 
summarized  by  Wolff  (1986).  Much  of  this  effort  has  been  directed  towards  the 
back-calculation  of  appropriate  values  for  the  analysis  parameters,  in  order 
to  verify  the  adequacy  of  the  levees  as  well  as  the  adequacy  of  design  crite¬ 
ria  for  future  levees.  To  back-calculate  any  parameter,  the  remaining  parame¬ 
ters  must  be  known  or  assumed;  back-calculated  parameters  can  only  be  as 
accurate  as  the  assumptions  made  for  the  "constrained"  parameters.  For  condi¬ 
tions  that  fit  the  assumptions  of  the  conventional  analysis,  the  relationship 
between  piezometric  elevation  at  the  levee  toe  and  the  river  stage  will  be 
linear  and  of  the  form 

y  =  mx  +  b  (2) 


or 

Piez.  el  =  (hQ  /H) (river  stage  -  ground  el)  +  ground  el  (3) 

The  slope  of  the  relationship,  iWH  ,  is  a  function  of  the  foundation  geome¬ 
try  (lengths  and  thicknesses)  and  permeability  ratio  (k^/K^).  The  intercept 
depends  on  the  ground  elevation  (or  elevation  of  landside  water) .  Where  the 
geometry  and  ground  elevations  assumed  in  the  analysis  match  those  in  the 
field,  then  the  field  permeability  ratio  can  be  determined  by  trial  and  error; 
the  correct  ratio  will  yield  a  straight  line  plot  through  the  observed  data. 
Where  the  real  ground  profile  is  complex,  assessment  of  the  field  permeability 
ratio  is  uncertain  because  of  the  uncertainty  regarding  the  appropriate 
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"effective  values"  of  the  layer  thicknesses  and  ground  elevations.  Different 
permeability  ratios  will  be  back-calculated  for  different  (although  reason¬ 
able)  assumptions  of  effective  layer  thicknesses  and  elevations. 

10.  By  using  the  programs  developed  herein,  the  foundation  geometry  is 
more  completely  described,  and  fewer  assumptions  are  required  to  make  the 
problem  tractable.  Thus,  back-calculated  permeability  ratios  should  be  more 
reliable  and  less  uncertain  than  those  obtained  from  conventional  analysis. 
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PART  II:  SPECIAL  CASES  OF  FOUNDATION  CONDITIONS 


11.  As  previously  stated,  the  conventional  assumptions  of  2-D  flow  and 
a  foundation  profile  of  two  uniform  horizontal  layers  are  not  always  consis¬ 
tent  with  actual  foundation  conditions.  For  this  research,  three  alternate 
sets  of  assumptions  have  been  identified  and  a  computer  program  for  the  finite 
difference  solution  of  steady  state  laminar  flow  in  porous  incompressible 
saturated  media  was  developed  for  each.  These  special  cases  are  representa¬ 
tive  of  many  prototype  locations.  In  many  cases,  the  relevant  deficiencies  of 
the  conventional  method  can  be  circumvented  by  selecting  an  appropriate 
alternative.  These  alternate  of  "special"  conditions  are  described  in  detail 
below. 

Foundations  Characterized  by  Three  Layers 

12.  The  special  condition  considered  is  the  case  of  a  foundation  con¬ 
sisting  of  three  materials  deposited  in  horizontal  layers.  The  pervious  mate¬ 
rials  immediately  below  a  clay  top  blanket  are  often  fine  sand  or  silty  sand, 
while  those  at  greater  depth  are  typically  medium  or  coarse  clean  sand.  Such 
layering  is  consistent  with  the  geologic  environment  of  a  meandering  stream, 
where  the  expected  grain  size  distribution  is  fining  upward.  In  the  developed 
analysis,  each  of  three  layers  may  have  different  thicknesses  and  different 
horizontal  and  vertical  permeabilities.  The  layers  are  herein  referred  to  as 
the  semipervious  top  blanket,  moderately  pervious  middle  stratum,  and  pervious 
substratum.  In  such  deposits,  the  top  blanket  may  be  a  backswamp  deposit,  the 
middle  stratum  may  be  a  point-bar  deposit  of  relatively  uniform  fine  sands  or 
silts,  and  the  pervious  substratum  is  representative  of  a  lower  section  in  the 
point  bar  or  a  "channel  lag"  deposit.  Figure  3  illustrates  a  three-layer 
foundation.  Use  of  the  conventional  analysis  requires  that  the  analyst  either 
convert  the  middle  stratum  to  a  relatively  thin  equivalent  layer  of  top  blan¬ 
ket,  or  consider  it  part  of  the  substratum  and  average  its  relatively  low 
permeability  with  higher  permeability  values  farther  down  the  profile. 

Jeither  of  these  approaches  may  adequately  model  flow  conditions  along  the 
base  of  the  top  blanket.  If  the  horizontal  permeability  of  the  middle  stratum 
exceeds  the  vertical  permeability  of  the  substratum,  residual  heads  at  the 
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Figure  3.  Example  of  a  three-layer  foundation 


base  of  the  top  blanket  may  be  quite  dependent  on  the  properties  of  the  middle 
stratum,  and  less  dependent  on  properties  of  the  substratum. 


Foundations  Characterized  by  Two  Layers  of  Irregular  Shape 


13.  The  top  blanket  thickness  and  the  ground  surface  elevation  are  not 
always  uniform.  In  most  cases,  prevailing  "average"  values  can  be  used  for 
analysis.  However,  in  certain  cases,  foundation  conditions  cannot  be  reason¬ 
ably  modeled  using  such  averages.  The  landside  ground  surface  may  slope  down¬ 
ward  away  from  the  levee;  this  is  commonly  the  case  where  flood-control  levees 
are  built  over  natural  levees.  Ridge  and  swale  topography  may  be  presented 
landward  of  the  levee.  Borrow  pits  and  ditches  interrupt  the  ground  surface 
and  form  significant  discontinuities  in  the  top  blanket.  The  occurrence  of 
sand  boils  is  highly  related  to  the  presence  of  discontinuities  such  as  those 
listed  above  (WES  1956a).  Where  clay-filled  channel  plug  deposits  run  paral¬ 
lel  the  levee,  the  top  blanket  may  be  locally  quite  thick,  and  seepage  will  be 
concentrated  at  the  edges  of  the  plug.  Some  of  the  factors  causing  irregular¬ 
ities  in  the  top  blanket  are  illustrated  in  Figure  4.  The  lower  Mississippi 
River  Valley  exhibits  markedly  irregular  top  blanket  conditions  because  of  the 
high  degree  of  stream  meandering  that  has  occurred.  Irregularities  in  the 
bedrock  surface  elevation  may  also  be  present,  due  to  features  such  as  pre¬ 
glacial  channels;  however,  seepage  conditions  are  usually  more  affected  by 
variations  In  top  stratum  thickness  than  variations  in  substratum  thickness. 


Angles  or  "Corners"  in  Levee  Alignment 


14.  Where  levees  parallel  major  river  bends  or  where  mainline  river¬ 
front  levees  meet  tributary  or  flank  levees,  the  levee  alignment  may  bend  and 
form  an  angle  or  corner.  In  such  cases,  the  area  of  the  levee  toe  may  be  sub¬ 
jected  to  concentrations  of  seepage  because  the  flow  regime  is  not  2-D.  Such 
conditions  are  illustrated  in  a  plan  view  in  Figure  5.  Although  such  condi¬ 
tions  can  be  found  throughout  the  Mississippi  River  Valley,  the  best  piezome¬ 
tric  data  for  such  conditions  were  found  in  the  St.  Louis  District.  In  other 
areas,  it  appears  that  straight  reaches  of  levee  have  been  deliberately 
selected  for  piezometer  locations  to  simplify  analysis. 
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Figure  4.  Example  of  irregular  foundation 


PART  III:  SELECTION  OF  LEVEE  REACHES  FOR  STUDY 


Selection  Criteria 


15.  Levee  reaches  in  the  Rock  Island,  St.  Louis,  Memphis,  and  Vicksburg 
Districts  were  screened  to  identify  reaches  where  actual  performance  could  be 
compared  to  predicted  performance  using  the  three  developed  computer  programs. 
Two  criteria  were  considered:  first,  the  reach  must  fit  the  particular 
foundation  conditions  of  interest  (three  layer,  irregular,  or  corner)  and 
second,  sufficient  and  reliable  performance  data  must  be  available.  The  most 
well-documented  performance  data  are  for  those  levee  reaches  reported  by  Cunny 
(1980)  for  the  Rock  Island  District,  those  reported  by  the  St.  Louis  District 
(USAED,  St.  Louis  1976),  those  reported  by  WES  (WES  1956a  and  b,  1964)  for  the 
Memphis  and  Vicksburg  Districts,  and  those  reported  by  McClelland  Engineers, 
Inc.  (1985)  for  the  Vicksburg  District.  Even  from  these  sources,  significant 
amounts  of  data  are  missing  or  have  been  identified  by  the  original  analyst  as 
having  questionable  reliability. 


Selected  Reaches 


16.  The  following  reaches  were  identified  in  the  preliminary  screening 
as  fitting  the  conditions  of  interest  and  having  reasonably  accurate  and 
complete  data.  Those  identified  by  an  asterisk  have  been  analyzed  for  this 
report  and  are  discussed  in  detail  in  Parts  IV  through  VI.  Locations  of  the 
reaches  analyzed  are  shown  in  Figure  6. 

Foundations  Characterized  by  Three  Layers 

*  Rock  Island  District,  Sny  Island  "F" 

Memphis  District,  Caruthersville 
Memphis  District,  Commerce 
Vicksburg  District,  Upper  Francis 

*  Vicksburg  District,  Eutaw 

Foundations  Characterized  by  Two  Layers  of  Irregular  Shape 
Rock  Island  District,  Sny  Island,  Range  "G" 

*  Rock  Island  District,  Hunt,  Range  "B" 

Rock  Island  District,  South  Quincy,  Range  "SQ" 

Rock  Island  District,  South  River,  Range  "SRC" 
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Figure  6.  Locations  of  analyzed  levee  reaches 


Foundations  Characterized  by  Two  Layers  of  Irregular  Shape  (Continued) 
St.  Louis  District,  Perry  County,  Sta  329+85 
Memphis  District,  Gammon 

*  Memphis  District,  Commerce 

*  Memphis  District,  Stovall 

*  Vicksburg  District,  Bolivar,  Range  "D" 

Angles  or  "Corners"  in  Levee  Alignment 

Rock  Island  District,  Bay  Island,  Range  C 
St.  Louis  District,  Columbia  Sta.  653 

*  St.  Louis  District,  Degognia  Sta.  260-290 
St.  Louis  District,  Grand  Tower  Sta.  430 

*  St.  Louis  District,  East  Cape  Sta.  94 
St.  Louis  District,  East  Cape  Sta.  309 
Memphis  District,  Farrell 

Memphis  District,  Stovall 


Vicksburg  District,  Bolivar,  Range  "D" 


PART  IV:  FOUNDATIONS  CHARACTERIZED  BY  THREE  LAYERS 


17.  To  analyze  underseepage  conditions  for  foundations  consisting  of 
three  layers,  a  computer  program  named  LEVEE3L  was  written.  Given  the  thick¬ 
nesses  and  effective  lengths  of  the  layers,  LEVEE3L  generates  a  grid  of 
96  points  or  nodes  (8  rows  by  12  columns)  and  solves  the  differential  equation 
for  2-D  flow  at  each  node  using  the  finite  difference  method.  The  spacing 
between  the  nodes  varies  in  both  the  x-  and  y-directions  and  is  a  function 
of  the  specified  geometry.  Figure  7  defines  the  variables  used  by  LEVEE3L  to 
describe  the  problem  geometry  and  illustrates  the  generated  grid.  The  ground 
surface  coincides  with  Row  1,  the  base  of  the  top  blanket  coincides  with 
Row  3,  the  base  of  the  middle  stratum  coincides  with  Row  5,  and  the  base  of 
the  substratum  coincides  with  Row  8.  Column  1  corresponds  to  an  open 
entrance,  and  Column  12  corresponds  to  an  open  exit.  Infinite  or 

distances  are  modeled  by  specifying  very  large  values,  as  is  traditionally 
done  in  finite  element  modeling.  The  landside  levee  toe  is  at  Column  7.  The 
exit  gradient  is  obtained  by  dividing  the  excess  head  at  node  (3,7)  by  the 
thickness  of  the  top  blanket.  The  finite  difference  equation  for  flow  at  an 
interior  node  is  shown  in  Figure  8.  Use  of  the  program  LEVEE3L  is  described 
in  Appendix  A.  Results  of  the  program  are  compared  to  the  other  programs  in 
Appendix  D. 


Effect  of  Moderately  Pervious  Middle  Stratum 


18.  The  effect  of  a  middle  stratum  was  initially  investigated  by  using 
the  program  to  perform  a  series  of  parametric  studies.  Results  of  these 
studies  are  illustrated  in  Figures  9  through  11.  For  all  parametric  studies, 
the  base  width  of  the  levee  was  arbitrarily  taken  as  50  ft  and  the  differen¬ 
tial  head  was  taken  as  20  ft.  Most  real  levees  of  20-ft  height  would  have  a 
wider  base  width;  this  would  tend  to  reduce  the  calculated  gradients  somewhat. 
In  the  first  study,  constant  foundation  geometry  was  assumed,  and  the  ratio  of 
the  vertical  permeability  in  the  top  blanket  (k^  )  to  the  horizontal  permea¬ 
bility  in  the  substratum  (k^)  was  assume^  to  be  1,000.  Then  the  exit  gradi¬ 
ent  through  the  top  blanket,  i  ,  was  evaluated  et  point  1 ,  7  in  Figure  7 
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Gradient 


Variables  and  grid  used  by  program  LEVEE3L 


varies 


permeability  of  middle  stratum 


v  Z1  “  10  *1V  ~  0*0002  ft/min  k..  varies 

Cl  Z2  varies  k2y  =  0.0010  ft/min  ^  varies 

0.70  -  d  “  6o'  k,v  varies  k3h  =  0.2000  ft/min 
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Figure  10.  Gradient  versus  thickness  of  middle  stratum 
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THICKNESS  OF  TOP  BLANKET,  z 


as  a  function  of  the  permeability  of  the  middle  stratum  and  the  permeability 
ratio  within  each  stratum.  The  results  are  presented  in  Figure  9.  It  is 
shown  that  the  gradient  increases  as  permeability  of  the  middle  stratum 
increases  and  the  gradient  decreases  as  the  permeability  ratio  within  the 
layers  increases. 

19.  In  the  second  study,  the  thickness  of  the  top  blanket,  thickness  of 
the  substratum,  and  the  vertical  permeabilities  were  fixed  and  the  gradient 
investigated  as  a  function  of  middle  stratum  thickness  and  permeability  ratio. 
The  results  are  presented  in  Figure  10.  It  is  shown  that  the  gradient 
decreases  as  the  middle  stratum  thickness  increases  and  the  gradient  increases 
as  the  permeability  ratio  within  the  layers  increases. 

20.  The  third  parametric  study  was  similar  to  the  second  except  that 
the  thickness  of  the  middle  stratum  was  fixed  and  the  thickness  of  the  top 
blanket  was  varied.  The  results  are  shown  in  Figure  11.  It  is  shown  that  the 
gradient  decreases  as  the  top  blanket  thickness  increases  (as  is  well  known) 
and  that  the  gradient  decreases  as  permeability  ratio  increases. 

Actual  Versus  Predicted  Performance 


21.  The  program  LEVEE3L  was  used  to  analyze  two  locations  where  a 
moderately  pervious  middle  stratum  is  present  and  where  piezometric  data  are 
available.  Although  LEVEE3L  is  quite  useful  to  predict  performance  where 
foundation  conditions  are  reasonably  well  known,  its  generality  leads  to  com¬ 
plications  when  used  to  back-calculate  permeability  values  or  rations.  Values 
must  be  assumed  for  five  unknowns  in  order  to  back-calculate  the  value  of  the 
sixth;  a  unique  solution  for  permeability  ratios  cannot  be  obtained.  Also, 
predictions  can  be  very  sensitive  to  the  selected  ground  elevation,  which 
often  must  be  assumed.  To  illustrate  the  capabilities  of  the  program,  it  was 

assumed  that  the  ratio  k,  :k,,  :k„  :k„,  :k„  :k„,  was  of  the  form  l:4:k_  :4k_  :- 

Iv  Ih  2v  2h  3v  3h  2v  2v 

250:1,000.  This  corresponds  to  a  ratio  of  1,000  in  a  conventional  analysis. 
The  program  was  then  used  to  back-calculate  the  values  of  k£v  and  k£^  .  A 
discussion  of  these  analyses  follows. 

Rock  Island  District,  Sny  Island  Range  F 

22.  This  piezometer  range  was  established  in  1954  on  the  east  bank  of 
the  Mississippi  River  at  river  mile  300.1  above  the  confluence  with  the  Ohio. 
It  was  previously  analyzed  by  Cunny  (1980).  A  soil  profile  for  the  site  is 
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shown  in  Figure  12.  The  top  stratum  thickness  ranges  from  4.8  to  10.0  ft  and 
generally  consists  of  2  to  4  ft  of  lean  clay  overlying  silt  and  silty  sand. 
Considerable  seepage  has  been  reported  at  this  location:  toe  seepage,  sand 
boils  and  pin  boils  in  1960,  pin  boils  in  1965,  and  light  toe  seepage  in  1973. 
These  observations  do  not  correlate  well  with  river  stage,  as  stages  were  the 
highest  in  1973  and  lowest  in  1960.  These  discrepancies  may  in  part  be 
related  to  levee  enlargement  between  1965  and  1967  and  provision  of  a  berm 
that  lengthened  the  effective  base  width. 

23.  Parameters  used  in  the  analyses  are  listed  in  Table  1.  The  column 
labeled  "conventional  analysis"  summarizes  values  presented  by  Cunny  (1980). 
Using  LEVEE3L,  the  profile  was  modeled  as  a  3.0-ft-thick  clay  top  blanket 
overlying  a  5.0-f t-thick  middle  stratum.  Two  computer  analyses  are  reported. 

In  analysis  "A,"  the  landside  distance  to  an  open  exit,  ,  was  taken  as 

3,000  ft  to  represent  a  foundation  with  infinite  landward  extent.  Analysis 
"B"  was  made  assuming  as  400  ft  to  match  the  observed  sluggish  response 

of  piezometer  F-4,  which  suggested  that  most  of  the  seepage  was  exiting 
relatively  close  to  the  levee.  Results  of  the  computer  analyses  are  plotted 
with  the  actual  data  in  Figures  13  through  16.  A  reasonable  match  to  the 
observed  performance  was  obtained  using  a  permeability  ratio  of  1:4:50: 200 : — 
250:1,000.  Although  the  finite  difference  grid  was  developed  primarily  with 
the  purpose  of  assessing  heads  and  gradients  at  the  levee  toe,  heads  and 
gradients  for  other  piezometer  locations  were  determined  by  interpolation  of 
the  final  heads  reported  in  the  program  output  file.  The  results  of  analyses 
"A"  and  "B"  are  virtually  identical  in  the  vicinity  of  the  levee,  but  analysis 
"B"  provides  a  better  match  remote  from  the  levee. 

24.  A  shortcoming  shared  by  b.  the  conventional  analysis  and  the 

three-layer  analysis  is  the  need  to  assign  a  single  value  for  the  ground  ele¬ 
vation  when  the  ground  is  in  fact  uneven.  The  appropriate  ground  elevation 
for  analysis  can  be  inferred  from  piezometric  response  using  Equation  3  in 
this  report,  as  the  residual  head  should  be  zero  for  a  river  stage  equal  to  the 
landside  ground  elevation.  A  value  of  457.5  was  determined  by  extrapolating  a 
line  through  the  piezometric  data  down  to  a  point  where  the  piezometric  eleva¬ 
tion  matched  the  river  elevation.  This  ground  elevation  is  2.0  ft  lower  than 
the  ground  elevation  at  piezometer  F-3  assumed  by  Cunny  (1980)  with  the  result 
that  higher  residual  heads  are  predicted  herein.  For  a  river  stage  of  472.8 
(levee  crest),  the  predicted  piezometric  elevation  at  F-3  is  460.8.  This 
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Table  1 

Parameters  Used  for  Analyses,  Rock  Island  District, 
Syn  Island,  Range  F 


Analysis 

Parameter* 

Conventional 
Analysis 
(Cunny  1980) 

LEVEE3L 
Computer 
Analysis  "A" 

LEVEE3L 
Computer 
Analysis  "B 

Lj  (ft) 

— 

510 

510 

l2  (ft) 

— 

100 

100 

l3  (ft) 

Infinite 

3,000 

400 

zl  (ft) 

6.9  ft  at  F-3 

3.0 

3.0 

z2  (ft) 

— 

5.0 

5.0 

d  (ft) 

34 

34.0 

34.0 

kly  (ft/min) 

— 

0.0002 

0.0002 

klh 

— 

0.0008 

0.0008 

k2v 

--- 

0.0100 

0.0100 

k2h 

— 

0.0400 

0.0400 

k 

3v 

— 

0.0500 

0.0500 

k3h 

— 

0.2000 

0.2000 

kf/ku 

31 

1,000 

1,000 

Levee  crest 

472.8 

472.8 

472.8 

Ground  el 

459.5 

457.5 

457.5 

s  (ft) 

219 

— 

— 

x3  (ft) 

108 

— 

— 

h  at  H 

1.3  at  F-3 

5.0  at  F-3 

3.3  at  F-3 

o  max 


*  Defined  in  Figures  2  and  7. 


32 


Piezometer  F-2,  1960,1 


3 


Figure  14.  Results  of  analyses,  Rock  Island  District,  Range  F,  piezometer 


Piezometer  F-3,  1960,1961,1962,1965 
Computer  Analysis  "A"  (LEVEE3L),  L3  =  3000 
Computer  Analysis  "B"  (LEVEE3L),  Lo  =  400 
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corresponds  to  a  head  of  3.3  ft.  As  the  tip  of  F-3  is  located  below  the  base 
of  the  middle  stratum,  the  exit  gradient  at  F-3  depends  on  the  assumption  made 
for  the  effective  top  blanket  thickness. 

25.  As  F-3  is  some  distance  landward  of  the  levee  toe  and  below  the 
middle  stratum,  it  does  not  directly  provide  information  regarding  seepage 
conditions  in  the  top  blanket  at  the  levee  toe.  The  computer  solution  gives  a 
residual  head  of  5.60  ft  at  this  point  with  the  river  stage  at  the  levee 
crest.  Dividing  this  value  by  a  top  blanket  thickness  of  3.0  ft,  a  gradient 
of  1.87  is  obtained,  well  above  critical.  This  is  consistent  with  past 
observations  of  heavy  seepage  and  boiling  at  Sny  Island  and  indicated  that  the 
boils  are  likely  related  to  the  presence  of  relatively  thin  clay  top  blanket 
deposits  underlain  by  silty  middle  strata  that  provide  little  seepage 
resistance . 

Vicksburg  District, 

Eutaw,  Miss.,  Line  D 

26.  This  piezometer  range  is  located  on  the  east  bank  of  the  Mississi¬ 
ppi  River  one-half  mile  from  the  town  of  Eutaw,  Miss.  The  soil  profile  at  the 
site  is  shown  in  Figure  17.  A  moderately  pervious  middle  stratum  of  5  to 

40  ft  of  very  fine  sand,  silty  sand,  and  sandy  silt  lies  between  the  top  blan¬ 
ket  and  pervious  substratum.  These  deposits  lie  in  an  ancient  channel  of  the 
Mississippi  River  (WES  1956a).  Piezometer  readings  are  available  from  1961 
(WES  1964)  and  1973  (McClellan  1985). 

27.  Parameters  used  in  the  analysis  are  listed  in  Table  2.  The  column 
labeled  "conventional  analysis"  summarizes  values  obtained  from  TM  3-424  (WES 
1956a).  Using  LEVEE3L,  the  profile  was  modeled  as  a  9.0-ft-thick  clay  top 
blanket  overlying  a  9.0-ft-thick  middle  stratum.  These  values  were  judgment- 
ally  selected  based  on  inspection  of  the  profile  and  the  fact  that  piezometer 
tops  are  located  in  a  number  of  different  materials.  Results  of  computer 
analyses  are  compared  to  observed  performance  in  Figures  18  through  22. 
Although  LEVEE3L  was  developed  primarily  for  analyzing  the  levee  toe,  other 
piezometer  locations  were  analyzed  by  interpolating  final  heads  in  the  program 
output  file.  A  reasonable  fit  through  the  scattered  observed  performance  in 
the  vicinity  of  the  levee  toe  (piezometers  D-6  and  D-7)  was  obtained  using  a 
permeability  ratio  of  1:4:25:100:250:1,000.  This  also  fit  the  limited  data  at 
riverside  piezometers  P-1  and  D-2 .  At  piezometers  D-3  and  D-4 ,  under  the 
levee,  the  slope  of  the  curve  is  approximately  correct.  Landside  of  the  levee 
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Cross  section  of  Vicksburg  District,  Eutaw,  Miss.,  Line  D  (after  WES  1964) 


Table  2 

Parameters  Used  for  Analyses,  Vicksburg  District, 

Eutaw,  Miss,  Line  D 

Analysis 

Parameter* 

Conventional 

Analysis 

LEVEE3L 

Computer 

Analysis 

Lj  (ft) 

2,500 

2,500 

l2  (ft) 

— 

450 

l3  (ft) 

Infinite 

750 

zl  (ft) 

18.0 

9.0 

Z2 

— 

9.0 

d  (ft) 

70.0 

70.0 

k^v  (ft/min) 

— 

0.00022 

k 

0.00088 

lh 

k 

_ 

0.0055 

2v 

k2h 

— 

0.0220 

k. 

3v 

— 

0.0550 

k3h 

0.2200 

0.2200 

kf/kbl 

800 

— 

Levee  crest 

161.0 

161.0 

Ground  el 

135.0 

135.0 

s  (ft) 

1,500 

— 

x3  (ft) 

1,000 

— 

h  at  H 
o  max 

10.4 

8.0 

*  Defined  in  Figures  2  and  7. 
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Piezometer  D-3,  1961,1973 
Piezometer  D-4,  1961,1973 

Computer  Analysis  (LEVEE3L) 
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Figure  19.  Results  of  analyses,  Vicksburg  District,  Eutaw,  Miss,  Line 

piezometers  D-3  and  D-4 


0  Piezometer  D-5,  1961,  1973 
A  Piezometer  D-6,  1961,  1973 
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Figure  20.  Results  of  analyses,  Vicksburg  District,  Eutaw,  Miss.,  Line 

piezometers  D-5  through  D-7 
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Results  of  analyses,  Vicksburg  District,  Eutaw,  Miss.,  Line 
piezometers  D-8  through  D-10 


Piezometer  D-11,  1961,  1 973 
Piezometer  D-12,  1961 


Vicksburg  District,  Eutaw,  Miss.,  Line 
ers  D-11  and  D-12 


toe,  the  predictions  fall  on  both  sides  of  observed  data,  depending  on  the 
piezometer  location  analyzed. 

28.  The  great  variability  in  piezometric  response  at  this  site  points 
up  the  problems  inherent  in  analyzing  irregular  soil  profiles.  Despite  the 
ability  of  LEVEE3L  to  account  for  the  thick  silt  layer  below  the  clay  top 
blanket,  the  analysis  remains  uncertain  and  subjective  because  the  irregular 
ground  surface  elevation  and  top  blanket  thickness  must  be  treated  as  single 
"effective"  values.  It  appears  that  the  program  LEVEEIRR  described  in  Part  V 
of  this  report  offers  some  advantage  over  LEVEE3L  for  such  conditions,  and  it 
may  be  desirable  to  develop  a  program  that  combines  some  of  the  features  of 
both . 
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PART  V:  FOUNDATIONS  CHARACTERIZED  BY  TWO  LAYERS  OF 

IRREGULAR  SHAPE 

Numerical  Modeling  Technique 

29.  To  analyze  underseepage  conditions  for  foundations  consisting  of 
two  layers  of  nonuniform  thickness  with  non-horizontal  boundaries,  a  computer 
program  named  LEVEEIRR  was  written.  LEVEEIRR  solves  the  same  differential 
equation  for  the  same  assumptions  as  Bennett's  (1946)  solution;  however,  the 
soil  boundary  elevations  and  the  layer  thicknesses  a  and  d  need  not  be 
constant  but  may  vary  as  a  piecewise  linear  function  of  horizontal  distance 

x  .  Input  to  the  program  consists  of  the  same  parameters  used  for 
conventional  analysis,  with  the  addition  that  the  elevations  of  the  top  of  the 
blanket,  top  of  the  substratum,  and  bottom  of  the  substratum  are  specified  for 
a  selected  set  of  x-coordinates .  Between  each  specified  x-coordinate,  the 
program  generates  nine  additional  nodes  along  the  interface  between  the  toi> 
blanket  and  substratum.  The  head  at  each  node  is  calculated  by  iteratively 
solving  a  finite  difference  equation  for  the  flow  at  each  node.  The  technique 
used  in  LEVEEIRR  is  illustrated  in  Figure  23.  The  program  output  provides  the 
residual  head  and  gradient  at  every  node  point,  and  thus  allows  calculation  of 
the  expected  gradient  at  local  discontinuities  such  as  landside  ditches,  thin 
spots  in  the  top  blanket,  and  edges  of  clay  plugs.  The  use  of  LEVEEIRR  is 
described  with  examples  in  Appendix  B. 

Effect  of  Nonuniform  Layer  Thickness 

30.  To  illustrate  the  capabilities  of  LEVEEIRR,  the  common  problem  of 
excavating  a  ditch  adjacent  to  a  levee  was  investigated.  It  was  assumed  that 
the  top  blanket  had  a  uniform  thickness  of  20  ft  before  excavating  the  ditch. 
If  was  further  assumed  that  the  thickness  of  the  substratum  was  100  ft,  the 
net  head  on  the  levee  was  20  ft,  and  the  permeability  ratio  was  500.  The 
depth  of  the  ditch,  D  ,  and  the  distance  from  the  ditch  crown  to  the  levee 
toe,  L  ,  were  then  varied  and  the  gradient  at  the  bottom  of  the  ditch  plotted 
as  a  function  of  the  two  variables.  Results  of  the  analysis  are  shown  in 
Figure  24.  As  would  be  expected,  it  is  shown  that  the  gradient  increases  with 
increasing  ditch  depth  and  decreases  with  increasing  distance  from  the  levee. 
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de  ditch 


Actual  Versus  Predicted  Performance 


31.  The  program  LEVEEIRR  was  used  to  analyze  four  levee  reaches  having 
markedly  irregular  subsurface  conditions.  For  each  reach,  the  assumed  perme¬ 
ability  ratio  was  varied  to  find  the  ratio  that  best  corresponded  to  observed 
performance.  A  discussion  of  these  analyses  and  results  follows. 

Rock  Island  District,  Hunt,  Range  B 

32.  This  piezometer  range  was  establisned  in  1957  on  the  east  bank  of 
the  Mississippi  River  about  25  miles  upstream  from  Quincy,  Ill.,  in  the  pool 
area  of  Lock  and  Dam  20.  It  was  previously  analyzed  by  Cunny  (1980).  A 
foundation  profile  and  the  idealized  section  used  for  computer  modeling  are 
shown  in  Figure  '.5.  Irregularities  in  the  profile  include  a  landside  ground 
elevation  about  5  ft  higher  than  the  riverside  and  an  impervious  top  stratum 
varying  from  5.3  to  8.2  ft  thick.  There  are  three  piezometers  at  the  site: 

B-l  near  the  levee  crest,  B-2  near  the  levee  toe,  and  B-3  on  a  small  ridge 
about  350  ft  landward  of  the  levee  centerline. 

33.  Parameters  used  in  the  analyses  are  shown  in  Table  3.  Piezometric 
elevations  calculated  using  LEVEEIRR  are  compared  to  observed  piezometric 
elevations  in  Figure  26.  Two  computer  analyses  were  performed.  Analysis  "A" 
assumed  a  permeability  ratio  of  64,  the  same  as  the  value  used  by  Cunny. 
Although  this  value  provided  a  reasonable  match  to  the  observed  data  at  the 
three  piezometer  locations,  a  somewhat  better  fit  was  obtained  by  reducing  the 
top  blanket  permeability  until  the  permeability  ratio  was  20.  The  latter 
results  are  labeled  analysis  "B."  It  should  be  noted  that  LEVEEIRR  is  capable 
of  predicting  piezometric  elevations  even  where  they  are  below  ground  level, 
as  is  the  case  for  B-l  and  B-3,  so  long  as  there  are  artesian  (confined  seep¬ 
age)  conditions  in  the  pervious  substratum.  Assuming  a  flood  level  at  the  top 
of  the  levee  (501.5)  and  an  interior  water  elevation  of  487. C,  a  maximum 
gradient  of  0.73  is  predicted  to  occur  at  the  landside  levee  toe  near 
piezometer  B-2. 

Memphis  District,  Commerce,  Miss.,  Line  H 

34.  This  piezometer  range  was  established  in  1942  about  10  miles  north 
of  Tunica,  Miss.  Heavy  seepage  damage  was  reported  at  the  site  during  the 
1937  high  water.  The  levee  is  located  about  2,200  ft  from  the  Mississippi 
River  on  terrain  characterized  by  numerous  ridges,  ditches,  and  swales.  A 
foundation  profile  and  the  idealized  cross  section  used  for  analysis  are  shown 
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Table  3 

Parameters  Used  for  Analyses,  Rock  Island  District 

Hunt,  Range  B 


LEVEEIRR 

LEVEEIRR 

Analysis  Parameter 

Conventional 

Computer 

Computer 

(after  Cunny  1980) 

Analysis 

Analysis  "A" 

Analysis  "B" 

Lj_  (ft) 

— 

310 

310 

l2  (ft) 

137 

140 

140 

l3  (ft) 

Infinite 

2,000 

2,000 

d  (ft) 

112 

Varies 

Varies 

z  (ft) 

5.3  to  8.2 

-(Varies;  5.0  at 

B-2)  - 

(ft/min) 

— 

0.1280 

0.1280 

(ft/mln) 

— 

0.0020 

0.0064 

kf/kbi 

64 

64.0 

20.0 

s  (ft) 

459 

— 

— 

x3  (ft) 

227 

— 

— 

Levee  crest 

501.5 

501.5 

501.5 

Ground  el 

489.5  at  B-2 

Varies 

Varies 

H 

max 

12.0 

-(Varies;  501.5  - 

ground  el)- 

h  at  H  at  B-2 

o  max 

2.7 

4.9 

3.7 

i  at  B-2 

max 

0.43 

0.99 

0.73 
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in  Figure  27.  Seventeen  piezometers  were  used  for  analysis  at  this  section. 
Piezometric  data  were  obtained  in  1961  (WES  1964).  As  some  of  the  piezometers 
are  located  on  different  but  adjacent  levee  cross  sections,  some  variation  is 
observed  for  piezometers  that  are  about  the  same  distance  from  the  levee. 

35.  Parameters  used  for  analysis  are  shown  in  Table  4.  The  column 
titled  "Conventional  Analysis"  provides  values  reported  in  TM  3-424  (WES 
1956a)  for  comparison.  Piezometric  elevations  calculated  using  LEVEEIRR  are 
compared  to  observed  data  in  Figures  28  through  32.  The  predictions  are  based 
on  a  permeability  ratio  of  514,  a  value  that  was  calculated  using  the  two 
permeability  values  reported  in  TM  3-424.  (Permeability  ratios  and  permeabil¬ 
ity  values  reported  in  TM  3-424  are  not  exactly  consistent) .  It  should  be 
noted  that  reasonable  predictions  can  be  obtained  using  LEVEEIRR  for  widely 
different  distances  from  the  levee  toe. 

Memphis  District, 

Stovall,  Miss.,  Line  B 

36.  This  piezometer  range  was  established  in  1948  about  3.5  miles  west 
of  Stovall,  Miss.  Seepage  damage  occurred  at  the  site  during  the  1937  high 
water.  A  foundation  profile  and  the  idealized  cross  section  used  for  analysis 
are  shown  in  Figure  33.  Irregularities  in  the  profile  include  a  very  thin 
riverside  top  stratum  due  to  removal  of  borrow  materials,  a  thick  landside  top 
stratum,  and  a  silt  plug  below  the  levee  centerline.  Piezometric  data  were 
obtained  in  1961  (WES  1964). 

37.  Parameters  used  for  analysis  are  shown  in  Table  5.  The  column 
titled  "Conventional  Analysis"  provides  values  obtained  from  TM  3-424  for  com¬ 
parison.  Two  analyses  using  LEVEEIRR  are  reported.  In  analysis  "A"  the 
permeability  ratio  was  taken  as  432,  the  ratio  of  the  two  individual  values 
reported  in  TM  3-424.  For  analysis  "B"  the  permeability  ratio  was  taken  as 
1,000.  Piezometric  elevations  calculated  using  LEVEEIRR  are  compared  to 
observed  data  in  Figures  34  through  36.  it  is  noted  that  the  computer  solu¬ 
tions  plot  below  the  observed  data  for  piezometers  in  the  substratum  (E-15  and 
E-17)  by  several  feet.  Relatively  high  piezometric  response  has  been  observed 
at  these  piezometers.  The  profile  at  Stovall  is  an  extreme  case  of  different 
blanket  materials  and  thicknesses  on  opposite  sides  of  the  levee.  It  is 
likely  that  the  top  blanket  permeabilities  are  also  much  different  on  the  two 
sides.  Attempts  to  match  the  response  of  E-15  and  E-17  by  decreasing  the 
riverside  blanket  thickness,  increasing  the  permeability  ratio,  and  shortening 
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Figure  27.  Actual  and  modeled  cross  sections  of  Memphis  District,  Commerce,  Miss.,  Line 


Table  4 

Parameters  Used  for  Analyses,  Memphis  District 
Commerce,  Miss,  Line  H 


Analysis 

Parameter 

Conventional 

Analysis 
(WES  1956a) 

LEVEEIRR 

Computer 

Analysis 

Lj  (ft) 

— 

1,850 

l2  (ft) 

— 

214 

l3  (ft) 

— 

3,365 

d  (ft) 

165 

Varies 

z  (ft) 

7.0 

Varies 

kf  (ft/min) 

— 

0.18000 

k^  (ft/min) 

— 

0.00035 

kf/kbi 

580 

514 

s  (ft) 

1,350 

— 

x3(ft) 

880 

— 

Levee  crest 

220.2 

220.2 

Ground  el 

197.5 

197.5 

H  (ft) 

max 

22.7 

22.7 

h  at  H  (ft) 

o  max 

9.0 

8.7 

55 


Piezometer,  H-6,  1961 
Piezometer,  H-7,  1961 
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Figure  28.  Results  of  analyses,  Memphis  District,  Commerce,  Miss.,  Line 

piezometers  H-6  and  H-7 


Piezometer  H-9V  1961 
Piezometer  H-1 i,  1961 
Piezometer  H-12,  1961 
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Figure  29.  Results  of  analysis,  Memphis  Districts,  Commerce,  Miss.,  Line 

piezometers  H-9,  H-ll,  H-12 
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Figure  31.  Results  of  analyses,  Memphis  District,  Commerce,  Miss 
Line  H,  piezometers  H-16,  H-17,  and  H-19 
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piezometers  H-20  through  H-22 


DISTANCE  FROM  LEVEE  IN  FEET 


6 1 


33.  Actual  and  modeled  cross  sections  of  Memphis  District,  Stovall,  Miss.,  Line 


Table  5 


Parameters  Used  for  Analyses 

,  Memphis  District 

Stovall,  Miss. 

Line  B 

Conventional 

LEVEEIRR 

LEVEEIRR 

Analysis 

Analysis 

Computer 

Computer 

Parameter 

(WES  1956a) 

Analysis  "A" 

Analysis  "B1 

ll  (ft) 

— 

7,110 

7,110 

l2  (ft) 

— 

243 

243 

l3  (ft) 

— 

2,657 

2,657 

d  (ft) 

40 

Varies 

Varies 

z  (ft) 

15 

Varies 

Varies 

(ft /min) 

— 

0.1900 

0.3000 

(ft/min) 

— 

0.00044 

0.0003 

kf/kn 

600 

432 

1,000 

s  (ft) 

800 

— 

— 

x3  (ft) 

750 

— 

— 

Levee  crest 

194.3 

194.3 

194.3 

Ground  el 

164.5 

164.5 

164.5 

H 

29.8 

29.8 

29.8 

max 

h  at  H 

14.4 

8.3 

8.3 

o  max 
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the  entrance  distance  all  had  relatively  little  effect.  A  modification  of 
LEVEEIRR  to  allow  different  permeabilities  on  opposite  sides  of  the  levee  may 
afford  a  better  analysis  of  such  sections. 

Vicksburg  District, 

Bolivar,  Miss.,  Line  D 

38.  This  piezometer  range  is  located  along  the  east  bank  levee  of  the 
Mississippi  River  2  miles  northwest  of  Benoit,  Miss.  The  river  at  this  site 
is  about  8  miles  from  the  levee;  however,  Bolivar  Chute  lies  about  1,200  to 
1,500  ft  riverward  of  the  levee.  A  line  of  nine  piezometers,  D-l  through  D-9, 
run  perpendicular  to  the  levee  at  this  site.  A  soil  profile  at  the  piezometer 
line  is  shown  in  Figure  37.  Irregularities  in  the  profile  include  riverside 
borrow  pits,  landside  sublevees,  a  landside  ditch,  and  a  massive  clay-filled 
abandoned  channel  about  1,000  to  2,000  ft  landside  of  the  levee.  The  effect 
of  this  channel  is  to  concentrate  seepage  between  the  landside  levee  tow  and 
the  channel. 

39.  The  analyses  gave  more  weight  to  1973  flood  data,  where  available, 
than  1961  data.  Data  for  1961  were  obtained  during  a  falling  river,  and  an 
inverse  relationship  between  river  stage  and  piezometric  elevation  was  appar¬ 
ent.  Data  points  at  higher  river  stages  are  generally  1973  data,  and  those  at 
lower  stages  are  generally  1961  data.  Data  from  piezometers  D-l  and  D-7 
appears  unreliable  as  it  plots  considerably  lower  than  other  piezometric  data. 
This  may  be  due  to  the  piezometer  tops  being  founded  in  fine-grained  blanket 
materials  rather  than  in  the  pervious  substratum.  Parameters  used  in  the 
analyses  are  shown  in  Table  6.  Results  of  the  analyses  for  all  nine  piezom¬ 
eter  locations  are  shown  in  Figures  38  through  41.  A  permeability  ratio 
(k^/k^)  of  1,000  was  found  to  provide  the  best  fit  to  the  observed  data. 

These  results  are  labeled  "Analysis  A."  Performance  predictions  for  an 
assumed  permeability  ratio  of  100  are  also  shown  (Analysis  B) ,  and  it  is  seen 
that  the  difference  is  relatively  small,  particularly  at  relatively  large 
distances  landside  of  the  levee. 

40.  Assuming  a  flood  to  the  project  flow  line  of  el  166.4,  a  landside 
tailwater  of  140.5,  and  a  permeability  ratio  of  1,000,  a  gradient  at  the  levee 
toe  of  0.82  is  calculated.  Assuming  a  permeability  ratio  of  100  reduces  the 
gradient  to  0.63. 
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DISTANCE  FROM  LEVEE  IN  FEET 


Lgure  37.  Actual  and  modeled  cross  sections  of  Vicksburg  District 

Bolivar,  Miss.,  Line  D 


Table  6 

Parameters  Used  for  Analyses,  Vicksburg  District 
Bolivar,  Miss.,  Line  D 


Conventional 

LEVEEIRR 

LEVEEIRR 

Analysis 

Analysis 

Computer 

Computer 

Parameter 

(WES  1956a) 

Analysis  "A" 

Analysis  "B" 

Lj  (ft) 

— 

2,500 

2,500 

l2  (ft) 

— 

340 

340 

l3  (ft) 

— 

2,150 

2,150 

d  (ft) 

90 

Varies 

Varies 

z  (ft) 

7.0 

Varies 

Varies 

(ft/min) 

0.2400 

0.2400 

0.0024 

(ft/min) 

— 

0.00024 

0.0024 

kf/kU 

100  to  200 

1,000 

100 

S  (ft) 

500 

— 

— 

x3  (ft) 

350 

— 

— 

Flow  line 

167.2 

166.4 

166.4 

Tailwater 

141.0 

140.5 

140.5 

H 

max 

26.2 

25.9 

25.9 

h  at  H  at  B-2 

o  max 

10.8 

10.6 

8.2 

i 

max 

>1 

0.82 

0.63 

Piezometer  D- 3»  19^1 
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Figure  39.  Results  of  analyses,  Vicksburg  District,  Bolivar,  Miss. 

Line  D,  piezometers  D-3  and  D-4 


Piezometer  D-5*  1961,  1973 
Piezometer  D-6,  1961,  1973 


Figure  40.  Results  of  analyses,  Vicksburg  District,  Bolivar,  Miss.,  Line 

piezometers  D-5  and  D-6 


Piezometer  D-7,  1961,  1973 
Piezometer  D-8,  1961,  1973 
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RIVER  ELEVATION,  FT,  NGVD 


PART  VI:  ANGLES  OR  "CORNERS"  IN  LEVEE  ALIGNMENT 


Numerical  Modeling  Technique 

41.  To  analyze  underseepage  at  angles  or  corners  in  levee  alignment 
where  flow  conditions  are  not  2-D,  a  computer  program  named  LEVEECOR  was 
written.  Input  to  the  program  consists  of  the  same  parameters  used  for 
conventional  analysis  plus  the  angle  and  radius  of  the  same  parameters  used 
for  conventional  analysis  plus  the  angle  and  radius  of  the  levee  bend  and  the 
distance  to  the  river  on  either  side  of  the  bend,  these  variables  are 
illustrated  in  Figure  42.  Using  the  specified  dimensions,  LEVEECOR  generates 
an  irregularly  shaped  grid  of  node  points  in  a  horizontal  plane  as  shown  in 
Figure  43.  The  grid  represents  the  entire  thickness  of  the  pervious 
substratum,  the  differential  equation  for  steady-state  flow  is  expressed  in 
finite  difference  form  at  each  node,  and  a  solution  is  obtained  by  iteration. 
Flow  is  assumed  horizontal  in  the  substratum  and  vertical  in  the  top  blanket. 
At  each  node,  the  horizontal  flow  in  the  substratum  from  fo'-’'  adjacent  nodes 
plus  the  vertical  flow  through  the  top  blanket  (downward  on  the  riverside  and 
upward  on  the  land  side)  must  sum  to  zero.  The  technique  used  is  approximate 
the  horizontal  flow  components  are  not  exactly  perpendicular  in  the  vicinity 
of  the  corner.  Gradients  between  nodes  are  approximated  as  the  difference  in 
head  divided  by  the  distance  between  nodes;  flow  areas  between  nodes  are 
approximated  using  the  average  distances  between  nodes.  The  analysis 
technique  is  illustrated  in  Figure  44.  To  ensure  that  results  are  generally 
consistent  with  the  conventional  solution  at  points  remote  from  the  corner, 
the  head  at  the  landside  boundary  (Row  1)  is  automatically  matched  to  the 
conventional  solution  for  a  foundation  of  infinite  landward  extent,  the  use 
of  LEVEECOR  is  described  with  examples  in  Appendix  C. 

Effect  of  Levee  Curvature 


42.  To  assess  the  effect  of  levee  curvature,  a  parametric  study  was 
performed  using  LEVEECOR.  Typical  dimensions  and  permeability  values  were 
assumed,  and  the  residual  head  at  the  levee  toe  in  the  corner  (node  3,8)  was 
plotted  as  a  function  of  the  bend  angle,  0  ,  and  the  radius  of  curvature. 
Assumptions  and  results  are  shown  in  Figure  45.  It  is  shown  that  there  is  a 
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Figure  43.  Grid  generated  by  program  LEVEECOR 


Flow  at  an  Interior  Node 


£q  to  (5,6)  "  0  "  q(5,5)  +  q(6,6)  +  q( 5,7)  +  q(4,6)  +  qDown 


Figure  44.  Flow  at  an  interior  node,  program  LEVEECOR 
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Figure  45.  Results  of  parametric  study  usi  .  LEVbhCOR 
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tendency  for  the  residual  head  to  increase  with  increasing  bend  angle,  but  the 
increase  is  relatively  small,  less  than  0.5  ft  for  the  example  investigated. 
Effects  may  be  more  pronounced  for  shorter  entrance  distances,  different 
permeability  ratios,  etc;  however,  time  did  not  afford  a  complete  investiga¬ 
tion  of  every  parameter.  As  the  radius  was  increased  to  500  ft,  the  head 
increased;  the  head  decreased  after  the  radius  exceeded  this  value.  Such 
variation  may  be  related  to  the  fact  that  the  model  geometry  changes  as  the 
radius  changes.  Ideally,  all  the  curves  should  pass  through  the  point  labeled 
"hand  solution";  however,  discrepancies  arise  because  the  numerical  solutions 
for  different  radii  generate  different  grid  geometries.  Some  degradation  of 
results  is  apparent  for  relatively  large  and  small  radii  and  for  bend  angles 
greater  than  90  deg. 

43.  Although  the  study  seems  to  indicate  that  the  bend  angle  may  have 
only  a  slight  effect  on  the  residual  head  for  a  symmetric  problem,  the  program 
LEVEECOR  should  be  useful  for  analysis  of  asymmetric  problems  because  of  its 
ability  to  model  different  distances  to  the  river  on  either  side  of  the  bend. 
LEVEECOR  should  be  used  for  radii  greater  than  500  ft  since  the  model  failed 
to  produce  reasonable  trends  or  agree  with  the  closed  form  solution. 

Actual  Versus  Predicted  Performance 


44.  The  program  LEVEECOR  was  used  to  analyze  a  set  of  corner  reaches 
where  piezometric  data  were  available.  Performance  predictions  obtained  from 
the  program  were  compared  to  predictions  based  on  conventional  analysis  and  to 
actual  piezometer  readings.  A  discussion  of  these  analyses  and  results 
follows . 

St.  Louis  District, 

Degognia,  Station  260-290 

45.  This  piezometer  range  is  located  in  the  Degognia  Levee  and  Drainage 
District  adjacent  tc  a  near-right-angle  bend  in  the  Mississippi  River  called 
Liberty  Bend.  A  plan  of  the  levee  and  a  foundation  profile  are  shown  in  Fig¬ 
ure  46.  On  the  upstream  side  of  the  bend,  the  Mississippi  River  is  about 

650  ft  from  the  levee.  On  the  downstream  side  of  the  bend,  a  chute  separates 
the  floodplain  from  Wilkinson  Island,  and  the  Mississippi  River  is  on  the 
opposite  sided  of  the  island.  I'nderseepage  performance  at  this  location  dur¬ 
ing  the  1973  flood  was  described  by  the  I'SAED,  St.  Louis  (  1  976). 
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ELEVATION,  FT,  NGVD 


Figure  46.  Flan  and  profile  oi  St.  Louis  District 
Degognia,  Station  d60-2‘>()  after  WHS  !V56a) 


300+00 


Piezometer  P-12  at  the  upstream  end  of  the  bend  flowed  over  a  4-ft  extension; 
the  estimated  head  exceeded  5.2  ft,  and  the  estimated  gradient  exceeded  0.43. 

46.  Three  analyses  were  performed,  one  using  the  conventional  method 
and  two  using  the  program  LEVEECOR:  analysis  assumptions  are  summarized  in 
Table  7.  The  conventional  analysis  is  based  on  parameters  obtained  and 
inferred  from  TM  3-430  (WES  1956b).  The  permeability  ratio  for  the  conven¬ 
tional  analysis  was  reduced  from  1,000  to  500  to  better  fit  the  observed  data. 
The  entrance  distance  of  1,000  ft  is  a  judgmental  average  of  the  widely  vary¬ 
ing  distances  on  either  side  of  the  bend.  Computer  analysis  "A"  was  performed 
using  the  same  permeability  ratio  as  the  conventional  analysis,  but  the  dis¬ 
tances  to  the  river,  L1A  and  LIB,  represent  prototype  conditions.  For  com¬ 
puter  analysis  "B"  the  permeability  ratio  was  varied  to  get  the  best  match 
with  the  observed  data.  Performance  predictions  from  the  three  analyses  are 
compared  to  observed  data  in  Figure  47.  For  the  dimensions  modeled,  it  was 
found  that  a  permeability  ratio  of  150  provided  the  best  match  to  observed 
performance.  Using  the  program  LEVEECOR  provides  a  convenient  method  for 
back-calculating  the  permeability  ratio  in  this  case  where  an  appropriate 
single  value  for  the  distance  to  the  river  required  for  conventional  analysis 
is  difficult  to  assess. 

St.  Louis  District, 

East  Cape,  Station  94 

47.  This  piezometer  range  is  located  at  an  abrupt  50-deg  bend  at  Sta¬ 
tion  94  of  the  East  Cape  Girardeau  levee,  which  is  across  the  river  from  Cape 
Girardeau,  Mo.  The  data  used  are  from  the  1973  flood  and  were  previously 
analyzed  by  USAED,  St.  Louis  (1976).  A  plan  view  and  foundation  profile  are 
shown  in  Figure  48. 

48.  Three  analyses  were  performed;  assumptions  are  summarized  in 
Table  8.  The  conventional  analysis  is  based  on  parameters  obtained  and 
inferred  from  TM  3-430  (WES  1956b).  Computer  analysis  "A"  was  performed  using 
the  same  permeability  ratio  as  the  conventional  analysis,  but  the  geometry  is 
more  completely  described.  For  computer  analysis  "B"  the  same  geometry  was 
used,  but  the  permeability  ratio  was  varied  to  get  the  beat  match  with  the 
observed  data.  Performance  predictions  from  the  three  analya.es  are  compared 
to  observed  data  in  Figure  49.  It  was  found  that  the  permeability  ratio  had 
to  be  lowerei  to  25  to  match  observed  performance.  The  USAED,  St.  Louis 
(1976;  reported  that  numerous  pin  boils  were  present  in  the  fields  landside  of 
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Table  7 


Parameters  Used  for  Analyses, 

St.  Louis  District, 

Degognia,  Station 

260-290 

Analysis 

Parameter 

Conventional 

Analysis 

LEVEECOR 

Computer 

Analysis  "A" 

LEVEECOR 
Computer 
Analysis  "B' 

lia  (ft) 

— 

5,000 

5,000 

lib  (ft) 

— 

650 

650 

s  (ft) 

1,000 

— 

— 

l2  (ft) 

— 

300 

300 

THETA  (deg) 

— 

110 

110 

RADIUS  (ft) 

— 

1,120 

1,120 

D  (ft) 

70 

70 

70 

z  (ft) 

14.0 

14.0 

14.0 

kf/kbi 

500 

500 

150 

x3  (ft) 

990 

— 

— 

Levee  crest  el 

382.7 

382.7 

382.7 

Ground  el 

360.7 

362.7 

362.7 

H 

22.0 

22.0 

22.0 

max 

h  at  H 
o  max 

9.06 

10.38 

7.8 

Table  8 


Parameters  Used  for  Analyses,  St.  Louis  District, 
East  Cape,  Station  94 


Analysis 

Parameter 

Conventional 

Analysis 

LEVEECOR 
Computer 
Analysis  "A" 

LEVEECOR 
Computer 
Analysis  "B' 

lia  (ft) 

— 

850 

850 

lib  (ft) 

— 

400 

400 

l2  (ft) 

— 

290 

290 

s 

850 

— 

— 

THETA  (deg) 

— 

50 

50 

RADIUS  (ft) 

— 

290 

290 

d  (ft) 

70 

70 

70 

z  (ft) 

12.0 

12.0 

12.0 

f — i 

as 

1,000 

1,000 

25 

x3  (ft) 

916 

— 

— 

Levee  crest  el 

358.5 

358.5 

358.5 

Ground  el 

340.3 

340.3 

340.3 

H 

18.2 

18.2 

18.2 

max 

h  at  H 
o  max 

9.5 

9.2 

4.0 

Piezometer  P-11 
Conventional  Analysis 
LEVEECOR  Analysis  "A" 


Results  of  analyses,  St.  Louis  District,  East  Cape,  Station  94 


this  location.  These  observations  are  consistent  with  the  unusually  low 
permeability  ratio  indicated. 
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PART  VII:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


49.  The  "conventional  method"  of  underseepage  analysis  is  restricted  to 
the  special  case  of  foundations  consisting  of  two  layers,  each  of  uniform 
thickness  and  with  horizontal  boundaries.  For  this  special  case,  there  is  a 
closed-form  solution  for  the  governing  differential  flow  equation.  Where 
actual  foundation  conditions  are  complex,  considerable  engineering  judgment 
must  be  used  to  assign  equivalent  or  "effective"  values  for  the  parameters 
describing  the  foundation  geometry.  Analysis  results  can  be  very  dependent  on 
the  judgments  made  and  may  vary  considerably  among  different  engineers  analyz¬ 
ing  the  same  soil  profile.  Using  numerical  methods  and  microcomputers,  the 
differential  equation  can  be  solved  for  more  general  cases  of  foundation  and 
levee  geometry.  This  reduces  uncertainty  introduced  in  the  analysis  regarding 
the  geometry  and  allows  the  analyst  to  focus  on  the  effects  of  other  vari¬ 
ables,  in  particular  the  material  permeabilities  and  landside  water 
elevations. 

50.  Three  microcomputer  programs  were  developed  to  perform  underseepage 
analysis  for  three  special  cases  of  foundation  conditions  that  are  representa¬ 
tive  of  many  locations  in  the  Mississippi  Valley.  Using  these  programs,  pre¬ 
dictions  of  piezometric  head  were  made  and  compared  to  observed  piezometric 
data  at  eight  locations.  The  field  permeability  ratios  were  estimated  by 
varying  values  of  the  program  input  parameters  until  predictions  matched 
observed  performance.  Although  the  programs  were  evaluated  by  reanalyzing 
past  performance,  they  would  be  equally  useful  for  new  analyses.  Because  grid 
size  effects  on  the  numerical  solution  were  not  parametrically  investigated 

.  the  grid  size  cannot  be  varied  by  the  user,  the  reader  is  cautioned  that 
exit  gradients  computed  by  these  programs  may  not  be  conservative.  Experi¬ 
ences  with  the  programs  ana  findings  from  these  studies  are  summarized  below. 

Foundations  charac- 
terized  by  three  layers 

51.  Where  foundation  sands  are  significantly  different  in  the  upper  and 
lower  parts  of  the  substratum,  a  three-layer  model  provides  a  more  consistent 
or  general  representation  of  prevailing  conditions  than  the  conventional  two- 
layer  analysis.  The  program  LEVEE3L  allows  modeling  of  such  conditions  and 
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allows  for  different  vertical  and  horizontal  permeabilities  in  each  of  the 
three  layers.  Some  of  the  effects  of  variations  in  the  middle  stratum  have 
been  demonstrated  by  parametric  studies. 

52.  The  program  LEVEE3L  should  be  useful  to  predict  underseepage  condi¬ 
tions  where  three-layer  conditions  are  known  to  exist  and  reasonable  estimates 
can  be  made  of  the  six  required  permeabilities.  A  second  case  of  a  three- 
layer  foundation  not  studied  herein  is  the  case  of  a  sandy  or  silty  top 
blanket  overlying  a  clay  layer.  LEVEE3L  would  be  applicable  to  such  condi¬ 
tions.  Because  of  its  ability  to  account  for  six  permeabilities,  LEVEE3L  is 
less  useful  for  back-calculation  of  permeability  values  or  ratios  because  five 
values  must  be  assumed  to  calculate  the  sixth. 

53.  The  analysis  at  Sny  Island  illustrates  how  a  thin  top  blanket  of 
clay  overlying  a  silty  middle  stratum  can  lead  to  boils  while  piezometers  at 
the  base  of  the  middle  stratum  indicate  only  moderate  residual  heads. 

54.  The  analysis  at  Eutaw,  Miss.,  illustrates  that  irregularities  in 
the  elevations  and  thicknesses  of  the  soil  deposits  may  dominate  the  problem 
more  than  the  effect  of  the  middle  layer.  Assigning  appropriate  effective 
ground  elevations  may  be  critical  to  accurately  predict  landside  heads,  and 
these  elevations  may  differ  significantly  at  different  points  along  the 
profile. 

Foundations  characterized 
by  two  layers  of  irregular  shape 

55.  Using  conventional  analysis  techniques,  the  calculated  residual 
head  is  a  function  of  both  the  assumed  equivalent  geometry  and  the  perme¬ 
ability  ratio.  Calculated  residual  heads  may  be  highly  uncertain  because  of 
the  uncertainty  introduced  in  the  representation  of  foundation  geometry. 

Using  the  program  LEVEEIRR,  a  two  layer  foundation  geometry  can  be  represented 
as  precisely  as  it  is  known.  Thus,  the  permeability  values  and  ratio  can  be 
isolated  as  che  independent  variables.  The  program  is  applicable  to  many 
common  conditions  of  irregular  foundation  geometry,  including  borrow  pits, 
ditches,  and  clay-filled  channels.  It  should  be  particularly  useful  in 
assessing  the  impact  of  proposed  changes  adjacent  to  levees,  such  as  construc¬ 
tion  of  new  ditches  or  elargement  of  borrow  areas,  problems  that  commonly 
occur  and  are  not  conveniently  assessed  using  conventional  methods.  A  partic¬ 
ularly  useful  feature  of  LEVEEIRR  is  its  capability  to  predict  or  assess  per¬ 
formance  (residual  head  and  gradient)  at  any  desired  distance  from  the  levee 
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or  any  piezometer  location  by  appropriate  specification  of  segment  boundaries. 
In  the  author's  opinion,  LEVEEIRR  probably  has  the  greatest  potential  of  the 
three  programs  for  extensive  use  in  routine  practice. 

56.  The  above  comments  on  calculating  residual  heads  apply  as  well  to 
back-calculation  of  permeability  ratios  from  observed  residual  heads.  As  the 
foundation  geometry  can  be  precisely  described  using  LEVEEIRR,  back-calculated 
permeability  values  should  be  more  reliable  than  those  related  to  "trans¬ 
formed"  geometry. 

Angles  or  "corners"  in  levee  alignment 

57.  Using  program  LEVEECOR,  underseepage  conditions  at  bends  in  levee 
alignment  or  levee  corners  can  be  analyzed.  Because  of  the  problem  complex¬ 
ity,  the  program  employs  more  approximations  than  the  other  two  programs.  For 
typical  problems  investigated,  the  increase  in  residual  head  due  to  the  corner 
effect  is  relatively  small,  typically  less  than  a  foot.  This  effect  alone 
would  not  seem  to  warrant  using  the  program  in  lieu  of  conventional  methods. 
However,  the  program  may  be  useful  in  practice  because  it  allows  the  analyst 
to  model  actual  distances  to  the  river  on  either  side  of  levee  rather  than 
assuming  an  "effective  distance"  that  somehow  averages  widely  different  condi¬ 
tions.  LEVEECOR  should  be  used  with  caution,  particularly  for  radii  less  than 
500  ft  since  the  model  failed  to  produce  reasonable  trends  or  agree  with  the 
closed  form  solution  in  that  range. 

Permeability  ratios 

58.  Permeability  ratios  cannot  be  reliably  back-calculated  using 
LEVEE3L.  Back-calculated  permeability  ratios  obtained  using  LEVEEIRR  and 
LEVEECOR  varied  from  25  to  1,000  and  were  typically  500  or  less.  The  low 
ratio  of  25  at  East  Cape  was  for  a  site  where  pin  boils  were  occurring,  indi¬ 
cating  that  the  apparently  high  top  blanket  permeability  was  due  to  open  path¬ 
ways  or  defects  in  the  blanket.  For  the  most  part,  the  ratios  calculated 
using  the  programs  were  consistent  with  those  obtained  by  experienced  analysts 
using  conventional  analysis.  However,  the  programs  allow  a  less  experienced 
analyst  to  model  the  geometry  and  foundation  conditions  as  they  are  perceived 
without  the  assignment  of  effective  dimensions. 

Comparison  of  programs  and  results 

59.  The  three  programs  and  the  conventional  method  were  used  to  analyze 
a  common  problem,  and  the  results  are  compared  in  Appendix  D.  Variation  among 
the  solutions  for  the  problem  analyzed  was  several  feet,  a  range  somewhat 
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larger  than  expected  at  the  outset  of  the  research.  Some  reasons  for  the 
variation  are  suggested  in  the  appendix,  they  include  the  use  of  fundamentally 
different  assumptions  in  deriving  the  flow  equations  and  truncation  error.  It 
is  expected  that  some  of  this  variation  could  be  reduced  by  using  finer  grids 
or  higher  order  approximations  for  the  derivatives.  Such  improvements  would 
require  a  second  iteration  in  program  development.  The  intent  of  the  present 
research  was  to  demonstrate  the  applicability  of  special-purpose  microcomputer 
programs  (as  compared  with  general-purpose  computer  programs)  for  the  problems 
of  interest.  Based  on  the  parametric  studies  performed  and  the  case  histories 
analyzed,  it  is  be1xeved  that  the  present  versions  of  the  programs  are  suit¬ 
able  and  should  be  recommended  for  comparison  studies  to  evaluate  the  effects 
of  borrow  areas,  ditches,  varying  permeabilities,  etc. 

Recommendations 


60.  Based  in  the  results  of  this  research,  the  following  recommenda¬ 
tions  are  made: 

a.  The  developed  programs  should  be  used  as  a  supplement  to  con¬ 
ventional  analysis  where  geometry  warrants. 

b.  Modifications  should  be  made  to  LEVEE3L  to  provide  for  variable 
grid  spacing. 

c.  Assuming  that  LEVEEIRR  would  be  the  most  used  of  the  three 
programs,  improvements  should  be  made  to  allow  for  different 
landside  and  riverside  blanket  permeabilities,  to  provide  for 
expedient  data  input  and  file  building,  to  provide  a  graphic 
display  of  foundation  geometry,  and  to  provide  for  more  than 
one  landside  water  level. 

d.  A  three-layer  version  of  LEVEEIRR  should  be  developed  that 
incorporates  features  of  LEVEE3L. 

e.  By  analyzing  additional  reaches  and  systematically  assessing 
all  results,  better  guidelines  should  be  developed  for  estimat¬ 
ing  the  field  permeability  ratio  as  a  function  of  material 
types  and  thicknesses. 

f.  The  extension  of  the  programs  to  analyze  reaches  where  relief 
wells  are  present  should  be  investigated. 

g.  The  programs  should  be  field  tested  to  determine  the  need  for 
additional  changes  and  improvements. 
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APPENDIX  A:  COMPUTER  PROGRAM  LEVEE3L:  UNDERSEEPAGE  ANALYSIS 
FOR  FOUNDATIONS  CHARACTERIZED  BY  THREE  LAYERS 


1.  The  program  LEVEE3L  performs  underseepage  analysis  for  levee  founda¬ 
tions  consisting  of  three  layers  of  uniform  thickness  with  horizontal  bound¬ 
aries.  The  program  was  written  in  FORTRAN77  and  runs  on  IBM  PCs  or  compatible 
microcomputers  using  the  MS-DOS  operating  system. 

2.  Input  to  the  program  is  from  a  seven  line  data  file  without  line 
numbers.  This  file  can  be  created  using  any  word  processing  or  text-editing 
program  that  produces  a  standard  ASCII  file.  An  example  input  file  is  shown 
in  Figure  A1 .  Program  variables  are  defined  in  Figure  A2. 

3.  The  program  is  executed  by  typing  the  command  LEVEE3L  with  the  file 
LEVEE3L.EXE  resident  on  the  default  drive.  The  program  will  then  ask  for  the 
name  of  the  input  file  and  the  name  of  the  output  file.  If  the  output  file 
already  exists,  it  will  be  written  over;  otherwise,  it  will  be  created.  A 
sample  run  of  LEVEE3L  is  shown  in  Figure  A3. 

4.  The  program  provides  output  to  two  devices:  the  console  (screen) 
and  the  output  file.  Output  to  the  screen  includes  the  numbe  of  iterations, 
the  residual  head  at  the  levee  toe,  and  the  gradient  at  the  levee  toe.  If  the 
user  desires  more  detailed  output,  the  output  file  can  be  displayed  or  printed 
using  any  word  processing  or  text-editing  program  that  works  with  standard 
ASCII  files.  The  output  file  contains  the  final  calculated  heads  at  each  node 
point  as  shown  in  Figure  A4. 

5.  The  program  source  listing  is  presented  in  Figure  5. 
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EXAMPLE  INPUT  FILE  FOR  LEVEE3L 


_ Values _ 

Example  Problem 

20.00  0.00 

1000.00  50.00  3000.00 

10. 0  20.0  60.0 
.0002  0.0100  0.0400 
.00100  0.0500  0.2000 
.005  500 


Variable  Names 
TITLE 
HR,  HL 
LI,  L2 ,  L3 
Zl,  Z2 ,  D 
KV1,  KV2,  KV3 
KHl ,  KH2,  KH3 
TOL  TRIES 


)les  are  defined  as  follows: 


descriptive  title  for  the  problem  (80  characters  maximum). 

distance  from  the  riverside  levee  toe  to  the  river  or  an  open 
ntrance. 

distance  from  the  riverside  levee  toe  to  the  landslde  levee  toe. 

distance  from  the  landslde  levee  toe  to  an  open  seepage  exit.  This 

very  long  (e.g.  twice  the  effective  exit  distance)  when  modeling 
ns  of  infinite  extent. 

elevation  of  water  on  the  riverside  of  the  levee.  If  the  landside 
landslde  ground  is  taken  as  the  datum,  HR  is  the  net  head  on  the 


elevation  of  water  on  the  landside  of  the  levee  or  the  elevation  of 
ground  if  no  water  is  present.  This  level  is  often  taken  as  the 
which  case  HL  is  0.0. 

thickness  of  the  top  blanket  in  feet  or  any  consistent  units. 

thickness  of  the  middle  stratum  In  feet  or  any  consistent  units. 
Figure  Al.  Example  input  file,  program  LEVEE3L  (Continued) 
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D  is  the  thickness  of  the  substratum  in  feet  or  any  consistent  units. 

KV1,  KV2,  and  KV3  are  the  vertical  permeabilities  of  the  top  blanket,  middle 
stratum,  and  substratum,  respectively,  in  feet  per  minute  or  any  consistent 
units. 

KHl,  KH2,  and  KH3  are  the  horizontal  permeabilities  of  the  top  blanket,  middle 
stratum,  and  substratum,  respectively,  in  feet  per  minute  or  any  consistent 
units . 

TOL  is  the  maximum  tolerance,  or  maximum  difference  in  head  at  any  node 
between  successive  iterations.  The  iteration  orocess  will  stop  when  the 
maximum  residual  is  smaller  than  TOL.  This  value  should  be  smaller  than  the 
desired  accuracy  of  the  solution  by  a  factor  of  10  or  more.  For  an  answer 
accurate  to  0.1  ft,  a  value  for  TOL  of  0.003  is  recommended. 

TRIES  is  the  maximum  number  of  iterations  that  the  program  will  be  allowed  to 
make.  It  is  provided  to  stop  the  program  in  the  event  of  non-converging 
solutions.  Generally  the  program  should  terminate  because  TOL  is  reached 
before  TRIES.  A  value  of  500  is  suggested. 


Figure  Al.  (Concluded) 
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Figure  A2.  Definition  of  variables,  program  LEVEE3L 


c> 

OLEVEE3L 

PROGRAM  LEVEE3L  - 

UNDERSEEPAGE  ANALYSIS  FOR  3 --LAYER  PROFILES 

Wrlttan  by  Thomas  F.  W olff  and  H.  A.  A1  -Moussawl 
Michigan  State  University 

For  the  U.  S.  Army  Corps  of  Engineers 
Release  1.0  September  1987 


XXXXXX 

X 

XXXXXX 

X 

XXXXXX 


XXXX  XXXXX  XXXX 

XX  XXX  XX 
XXXXXXXI ixxxxxx 

XXI  IXI  IXI  IX 
XXXXXXXI  IXXXXXX 


ENTER  INPUT  FILE  NAME 
DATA3L 

ENTER  OUTPUT  FILE  NAME 
0UT3L 


Ex amp 1 •  Pr  ob 1  am 


INPUT  DATA 

RIVERSIDE  HEAD  LANOSIDE  HEAD 

20, 00  . OO 


LI 

L2 

L3 

1000.00 

30.00 

3000.00 

Zl 

Z2 

D 

10.00 

20.00 

60.00 

KYI 

KY2 

KY3 

KX1 

KX2 

KX3 

RATIO 

RATIO 

RATIO 

.000200 

.010000* 

.040000 

.001000 

.050000 

. 200000 

3.00 

IX 

O 

O 

5.00 

TOLERANCE  MAX  ITPra^IONS 

. 00300  500 


JUST  A  MOMENT,  I  AM  THINKING 

SOLUTION  COMPLETE  AFTER  191  ITERATIONS 
MAXIMUM  RESIDUAL  -  .00498 


RESIDUAL  HEAD  AT  BASE  OF  TOP  BLANKET  -  9.6 U 

GRADIENT  THROUGH  TOP  BLANKET  -  .97 

HEAD  LOSS  THROUGH  MIDDLE  STRATUM  -  .35 

GRADIENT  THROUGH  MIDDLE  STRATUM  -  .02 


OUlPoi  Z>-  "~Z  IN  FILE  OUT 3L 


WANT  TO  RUN  ANOTHER  PROBLEM  7  Y  OR  N 
N 

Stop  ~  Program  tarminatad. 

C> 


Figure  A3.  Example  run  for  program  LEVEE3L 
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EXAMPLE  OUTPUT  FILE  FOR  LEVEE3L 


Example  Problem 


INPUT  DATA 

RIVERSIDE  HEAD 

20.00 

LI 

1000.00 

Z1 

10.00 

KYI 

XXI 

RATIO 

.000200 

.001000 

5.00 

TOLERANCE 

.00500 

INITIAL  HEADS 


20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

.00 

.00 

.00 

.00 

00 

20.00 

16.67 

15.00 

13.33 

11.67 

10.00 

0.33 

6.67 

5.00 

1.67 

.00 

20.00 

16.67 

15.00 

13.33 

11.67 

10.00 

8.33 

6.67 

5.00 

1.67 

.00 

20.00 

16.67 

15.00 

13.33 

11.67 

10.00 

6.33 

6.67 

5.00 

1.67 

.00 

20.00 

16.67 

15.00 

13.33 

11.67 

10.00 

8.33 

6.67 

5.00 

1.67 

.00 

20.00 

16.67 

15.00 

13.33 

11.67 

10.00 

8.33 

6.67 

5.00 

1.67 

.00 

20.00 

16.67 

15.00 

13.33 

11.67 

10.00 

8.33 

6.67 

5.00 

1.67 

.00 

20.00 

16.67 

15.00 

13.33 

11.67 

10.00 

8.33 

6.67 

5.00 

1.67  .00 

SOLUTION  COMPLETE  AFTER  191  ITERATIONS 
MAXIMUM  RESIDUAL  -  .00*98 

FINAL  HEADS 

20.00  20.00  20.00  20.00  20.00  20.00  .00  .00  .00 


LANDSIDE  HEAD 

.00 


L2 


50.00 


L3 

3000.00 


22 


20.00 


60.00 


JCY2 

KX2 

RATIO 


XV  3 
KX3 
RATIO 


.010000 

.050000 

5.00 


.040000 

.200000 

5.00 


MAX  ITERATIONS 
500 


Figure  AU.  Example  output  file,  program  LEVEE 3L 


.00 

3.33 

3.33 

3.33 

3.33 

3.33 

3.33 

3.33 


.00 

(Continued) 
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20.00 

19.24 

18.46 

17.62 

16.68 

15.44 

4.88 

2.65 

1.66 

.51 

,00 

20.00 

18.48 

16.92 

15.24 

13.36 

11.08 

6.68 

5.30 

3.32 

1.03 

.00 

20.00 

18.45 

16.86 

15.14 

13.23 

10.92 

9.86 

5.40 

3.39 

1.05 

.00 

20.00 

18.43 

16.80 

15.06 

13.10 

10.80 

10.03 

5.50 

3.45 

1.07 

.00 

20.00 

18.41 

16.78 

15.02 

13.05 

10.75 

10.10 

5.54 

3.48 

1.07 

.00 

20.00 

18.41 

16.76 

15.00 

13.03 

10.74 

10.15 

5.56 

3.49 

1.08 

.00 

20.00 

18.41 

16.76 

14.99 

13.01 

10.74 

10.17 

5.57 

3.50 

1.08  .00 

RESIDUAL  BEAD  AT  BASE  OF  TOP  BLANKET  -  S.E8 
GRADIENT  THROUGH  TOP  BLANKET  -  .97 

HEAD  LOSS  THROUGH  MIDDLE  STRATUM  -  .33 

GRADIENT  THROUGH  MIDDLE  STRATUM  -  .02 


Figure  A4.  (Concluded) 


1.04 

2.09 

2.13 

2.17 

2.18 

2.19 

2.19 
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PROGRAM  LEVEEJl 


cmtmfmmmiitftmimtfmmmmimtmmmmimm 


c  t 

C  THIS  PROGRAM  HAS  HRJTTEN  BVi-  ( 
C  DR.  THOMAS  F.  NOLFF)  ASSISTANT  PROFESSOR  $ 
C  AND  t 
C  HAaaAN  It.  ALNQUaaAHlj  PH.D  CANDIDATE  t 
C  ( 
C  USING  THE  A.  H.  CASE  CENTER  t 
C  COMPUTER  AIDED  ENGINEERS  FACILITY  AT  I 
C  MICHIGAN  STATE  UNIVERSITY  ,  1987  t 
C  $ 


c 

C  VERSION  i.O  SEPTEMBER  10,  1987 
C 

CHHIHIttHHHinHIIIHtMliHIHIHIIIIIKHIIttllHimiimm 

C 

C  THIS  PROGRAM  USES  THE  FINITE  DIFFERENCE  METHOD  TO  ANALYZE  SEEPA6E 
C  THROUGH  A  LEVEE  FOUNDATION  CONSISTING  OF  THREE  LAYERS,  EACH  HUH 
C  DIFFERENT  HORIZONTAL  AND  VERTICAL  PERMEABILITIES. 

C 

CtHHtHiHlitHHIHIHItlHmHHHmillHHHtHtHimtllHtll 

C 

C  VARIABLES 
C 

C  6R1  Gradient  at  Inn  toe  through  top  blanket 

C  GR2  Gradient  at  levee  toe  through  aiddie  itratua 

C  H(R,C>  Head  at  row  R  and  coluan  C 

C  HI  Head  an  landside  of  levee  (ground  or  ta) 

C  HR  Head  on  riveriide  of  lave* 

C  KK1  Horizontal  peraeability  In  top  blanket 

C  KI2  Horizontal  peraeability  in  aiddie  itratua 

C  KX3  Horizontal  peraeability  in  tubitratua 

C  KYI  Vertical  peraeability  in  top  blanket 

C  KY2  Vertical  peraeability  in  aiddie  itratua 

C  KY3  Vertical  peraeability  in  tubitratua 

C  Li  Distance  froa  RS  levee  toe  to  river 

C  L2  Base  width  of  levee 

C  L3  Distance  (roa  LS  levee  toe  to  point  of  zero  head 
C  MR  Maxiaua  residual 

C  PR1  Peraeability  ratio  in  top  blanket 

C  PR2  Peraeability  ratio  in  aiddie  itratua 

C  PR3  Peraeability  ratio  in  substratua 

C  TOL  Tolerance 

C  TRIES  Maxiaua  nuaber  of  iterations 

C  Zi  Thickness  of  top  blanket 

C  Z2  Thickness  of  aiddie  stratu* 

C  Z3  Thickness  of  substratua  (pj 


Figure  A5.  Listing  of  program  LEVEE3L  (Sheet  1  of  8) 
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CHARACTERS  F1LEIN 
CHARACTER *44  FILEOUT 
CHARACTER*80  TITLE 
CHARACTER*!  AMS 
REAL  LI, L2.L3.HR, HL,Z1,Z2,73 
REAL  Kn,K)t2,Kl(J,KYllKY2,m,T0L 
REAL  DX1 ,DX2,DX3,D21  fDZ2,DZ3 
REAL  H1,H2,I13,N4,H5,H6,NR,NRN 
REAL  H!8,!2),0LDHI8,12) 

REAL  PR1,PR2,PR3,RHD1,RHD2,8R1,GR2,CFL0AT 
INTEGER  R,C,TRIES,ITT,I,J 

C  DISPLAY  INTRODUCTION 
CALL  INTRO 

C  OPEN  DATA  FILES 

100  MRITE  <*,H  '  ENTER  INPUT  FILE  HAHE* 

READ  (*, ' (A) ')  F1LEIN 

OPEN (ONI T'10,FILErFILEIN,  ERR=100,  STATUS*'OLD') 
REMIND  (10) 

NRITE  (»,»)  1  ENTER  OUTPUT  FILE  NAME' 

READ  («,  '(A) 'I  FILEOUT 
OPEN  (UNIT320,F1LE=FILE0UT) 

REMIND  (20) 

C  READ  DATA 

REA5(I0,*(A)')TITLE 

READdO,*)  HR,HL,L1,L2,L3,Z1,Z2,I3,KY1,KY2,KY3, 
1KX1 ,KX2,KX3, TOL , TRIES 

C  INITIALIZE  HEADS 

DO  MO  0=1,6 
140  H ( 1 ,CT=HR 

DO  130  C*7,!2 
150  H(l,C)=HL 

DO  140  R=2,8 
H (R , l ) =HR 
160  H(R,12)-HL 

DO  170  C-2, 1 1 
DO  170  R=2,B 
CFLOAT=C 

170  H(R,C)=HR-!CFL0AT/12)»(HR-HL) 


C  THE  PERMEABILITY'  RATIOS  ARE 
PR1 =KX1 /KYI 

Figure  A5.  (Sheet  2  of  8) 


A9 


PR2=KX2/KY2 

PR3-KX3/KY3 


C  PRINT  INPUT  DATA  AND  WRITE  TO  FILE 

WRITE(20,«)  TITLE 
WRITE (*,«) 

WRITE(«,*)  TITLE 

WRI TE  (20 ,910) 

NRITE(«,9I0) 

910  FORHAT!/'  INPUT  DATA') 

WRITE (20,9t  2) 

WR1 TE (* ,912) 

912  FORMAT ( / ‘  RIVERSIDE  HEAD  LANDSIDE  HEAD'! 

WRITE (20 , 91 4)  HR,HL 
WRITE («,914>  HR ,HL 
914  FQRHAT(2(F9.2, 111)1 

WRITE (20,916) 

WRITE!*, 916) 

916  FORMAT!/'  LI  12  LJ’J 

WRITE (20,918)  ll,L2,LJ 
WRITE!*, 918)  L1,L2,L3 
918  FORMAT (1X,3(F9.2,11X)) 

WRITE (20,920) 

WRITE!*, 920) 

920  FORMAT!/'  II  22  D) 

WRITE ( 20,922)  Z1,I2,23 
WRITE («,922)  Zl,22,23 
922  FORMAT (1X,3 (F9.2, 1 IX) ) 

WRITE(20,924) 

WRITE!*, 924) 

924  FORMAT!/ '  KYI  KY2  KY3' 

1/'  KX1  KX2  KI3' 

21'  RATIO  RATIO  RATIO') 

WRITE (20,926)  KY1,KY2,KY3,KXI,KX2IKI3,PR1,PR2,PR3 
WRITE!*, 926)  KYl,KY2,KY3,Kn,KX2,K!!J,PRl)PR2lPRJ 
926  FORMAT  </ ,1X,3!F9.6,I1X)  ,/,lX,3(F9.6,llX)  ,/,lI,3(F9.2,UI)) 

WRITE (20,928) 

WRITE!*, 928) 

928  FORMAT!/'  TOLERANCE  MAI  ITERATIONS') 

NRITE(20,93QI  TOL, TRIES 
WRITE!*, 930)  TOL, TRIES 
930  FORMAT  <1X,F9.5,1SX,IS) 

WRITE (20,932) 

932  FORMAT!/'  INITIAL  HEADS') 

WR I TE 1 20 , 934 1  ( (H(R,C) ,C=1 ,12) ,R=1 ,8) 
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C  CALCULATE  DIFFERENTIAL  DISTANCES 

DZl=Zl/2 
022=12/2 
DI3=ZJ/J 
0XI=Ll/5 
DX2=L2 
DI3=L3/5 
N1=0XI/D2J 
M2=DX3/DZ3 
H3=DX1 /DZ2 
N4=DX3/DZ2 
H5=DX1/DZ1 
H6=DX3/DZ I 

C  KitttlHHdHHiMiKHHMtlH 
C  ITERATION 

C  mHHHIUtHlltmtimlHHH 

C  INITIALIZE 
ITT  *0 

WRITE!*, •) 

WRITE!*,*)  '  JUST  A  NONENT,  I  AN  TH1NKIN6' 

C  BE6IN  NEXT  ITERATION 

300  ITT=HT+1 

C  WRITE  <»,30l)  ITT 

C  301  FORMAT (IX’  ITERATION  15) 

IFUTT+1.6T, TRIES)  60  TO  2400 

C  RESET  OLD  HEADS 

DO  400  R=l,8 
DO  400  01,12 
400  OLDH(R,C)=HIR,C) 

C  ROM  "1BHT  -  IHPERHEBLE  BOUNDARY 

DO  500  C=2,5 

500  H(8,C)*!!(KX3*!H!8,C-l)*H!8,C*l)>>/Ml)+t2*KY3*H17,C)*Ml))/ 

1  I2*KX3/B1*2«KY3*MI) 

DO  600  C=8 , 1 1 

600  H<B,C)*<<(KX3t(H(8tC-lHH(8,01>)>/N2M!2»KY3»HI7,C>*H2))/ 

1  (<2«KU1/N2*2«KY3*H2) 

H!8,&)=!.5*KY3«(<DX2*DXi)/DZ3)«H(7,6)*KX3*!(DZ3/DXl)* 

I  H!8,3)*(DZ3/DX2)*H!8,7)1 )/!.5*KY3*l !DX2+0X1)/DZ3)*KX3* 

1  (DZ3/DX2+DZ3/DXD) 

H!8,7)=!.5*KY3*!I0X2*DX3)/DZ3)*H(717) *KX3«! (DZ3/DX2) ♦ 
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1  H(8,6)  +  (D23/DI3)tK(8,6)l)/(.5«KY3*((DX2*DX3)/DZ3)*M3* 

1  (DZ J/DX2^DZ3/0X3) I 

C  ROMS  SIX  AND  SEVEN  -  INTERIOR  NODES 

00  700  1=6,7 
00  700  C=2,5 

700  H<I,C>=U(KI3»(H(I,C*l)*H(l,C-l>))/Nl)*lKY3MH(M,C)* 

I  HU*l,CI)«Nl)l/<2«ia(3/Nl*2*KY3»«l) 

00  800  J=6,7 

00  800  C=8,il 

800  H<J ,C)  =  ( ( (KX3*(H(J,C+1)*H(J,C-1>  »)/H2)  +  (KV3t (H(J-1 ,0 ♦ 

1  H(J*1,C))*N2))/(2*KX3/H2*2*KY3*I12) 

DO  900  1=6,7 

H(I,6)«IKX3»(DZ3/0X2UH(I,7>*KX3«(0Z3/DII1HH(I,5)*KY3* 

1  ( (0X1*0X21 /(2»DZ3))*(H(l-i, 6) *H(I*1, 61  ))/(KI3«(0Z3/ 

1  DXl+DZ3/DX2>+KY3*( <DXl4DX2>/DZ3) ) 

H(I,7)=>(KX3*(DZ3/0X3)*H(r,8)+(CI3*(DZ3/DX2)»H(I,6!+KY3* 

I  ( 10X2*0X3) / (2«DZ3) >«(H( t-1 ,7)*H(I*1 ,7) ) ) / (KX3M0Z3/ 

1  DX3+DZ3/DX2)+KY3*(IDX2*DX3)/DZ3)) 

900  CONTINUE 

C  ROM  FIVE  -  INTERFACE  NODES 

DO  1000  C=2,5 

1000  H<5,C>*im*(OXl/OZ3)iHI6,C)*KY2*<DXl/OZ2>*H<4,CI* 

1  ( ( (KX3+KX2) /21*< (DZ2+DZ3) 72) *H (5,C+l) ) /DX1-* » IKX3+KX2) / 

1  2) *((DZ2*DZ3)/ 12*0X1) )»H(5,C-l>)/(KY3*(DXl/DZ3)*XY2* 

1  <DX1/DZ2)+(((KX3+KX2)/2)*(0Z2*DZ3))/DI1) 

DO  1100  C=8, 11 

1100  H(5,C)=(<KY3*IOX3/OZ3)*H(6,C)*KY2*<OX3/DZ2I*H(4,C)* 

1  (((KX3+KX2)/2)»((DZ2+DZ3)/2)*H(5,C+l))/DX3*(tKX3+KX2)/ 

1  (2*0X3n*H(5,C-l>«(DZ2*0Z3>/2)>/<<m*<0I3/DZ3))*KY2* 

1  (DX3/DZ2)+(t<KX3+KX2)/2)*(0Z2+DZ3>)/DI3) 

H(5,6)=(XY3*H(6,6)*((DX1*DX2)/DZ3)+KY2*((DX1+DX2)/ 

1  DZ2)«H(4,6)*(<KX3*KX2)/2)t.5»<(DZ2*DZ3)/DX2)*H(5,7>* 

1  ( (KX3+KX21 /2)*.5*( (DZ2+0Z3) /DX 1 1 *H (5 , 5) ) / ( (KY3/DZ31* 

1  (DX1>0X2) ♦ (KY2/DZ2) •(DX1 ♦DX2) ♦ ( ( (KX3+KX2 ) 721  *.5* 

1  (DZ2+DZ3)  1/DX2*(( (KX3*KX2172)*.5*(DZ2*DZ31 )/DXll 

H(5,7) *(KY3*H (6,7) *((DI2*DX3)/DZ3)*KY2«( (0*2*0X31/ 

1  0Z2)*HI4,7)*((KX3*KX2)/2)O*(  ID22*DZ3)/DX3) *H(5,8)* 

1  ((KX3*KX2)/2)*.5*((DZ2*DZ3)/DX2)*H(5,6)I/<(KY3/DZ3)» 

1  (DJ2*DX3)*(KY2/0Z2)*(0*2*0*3)*( ( (KX3*KX2)/*)».5* 

1  (DZ2*DZ3) )/DX3*( ( (KX3*KX2)/2)*.5*(DZ2*DZ3) 1/DX2) 

C  ROM  FOUR  -  INTERIOR  NODES 

DO  1200  C=2 , 5 
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1200  H(4)C)  =  ((KX2t(HlA,C*l)*HM1C-U)l/l1J*KY2MHl3,C)*H(5,Cl)« 

1  H3)  /  (2*KX2/H3*2«KY2*I13) 

DO  MOO  C-8,11 

MOO  HM,CI  =  (IKX2»(H(4,01HHM,C-l)))/H4+KY2*lH(J,C)  *H(5,C)  I* 

1  H4)/<2*KX2/H4*2«’aY2«IH) 

H(4,6)S(KX2*(DZ2/DX2)*HM,7MKX2*(DZ2/0X1)«H(4,5MKY2» 

1  ((DJ1+0X2)/(2»DZ2I)«(H(3IM*H(5,6I)I/(U2»(3Z2/DY2+DZ2/ 

1  DX1 )  ♦  (KY2*  (DXM0I2)  )/0Z2> 

H(4,7M(t!;2»(DZ2/0X3)»HM,8MKX2*(DZ2/D«2)»HM,4)^Y2* 

1  (<DX2*DX3}/(2*DZ2))tlHI3,7)*H(5,7)))/(KX2*ID22/DX3*DZ2/ 
l  DX2M(KY2«(DX2*DX3)I/DZ2) 


C  RON  THREE  -  INTERFACE  NODES 

DO  1500  C=2,5 

1500  H(3,C>«(KY2*(DX1/0Z2MHM,CHKYM(0X1/DZIMH(2,CM 
1  (((KX2YKXl)/2)»((DZMDZ2)/2MH(3,C+m/DXMUKX2+KXl)/ 

1  2M(((DZl+DZ2l/2)/DXll*H(31C-l)l/(KY2«(0n/D22Mt;YM 
1  DX1/DZH(((KX2*KXI)/2)»(DZMDZ2))/DX1) 

DO  MOO  C=8,ll 

MOO  HI3,C)*(KY2»(DX3/DZ2)»HI4,C)+KYM(DX3/DZ1)*H(2,CM 
1  (UNX2+KX1)/2)*<:DZMDZ2)/2)*H(3,CH))/DX3M(KXMI(X2)/ 

1  2)*<<DZM022)/(2»DX3)MH(3,C-in/(KY2*(0X3/DZ2HKYM 
1  (DXJ/0Z1  W  <  mm  2)/2)HDUm2) )  /DU) 

H(3,6)MKY2»H(4,6MUDXMDX2)/DZ2MXYM((DIMDX21/ 

1  DZlMH(2,6)  +  ((KXl+tCX2l/2)«.5*<<DZMDZ2)/DX2MH(3,7M- 
1  ((KX2+KX1)/2)«.5KIDZMDZ2)/DX1MHI3,5))/((KY2/DZ2)« 

1  (DXMDX2MKYMt(0XUDI2)/DZlM((IUMKX2)/2)*.5» 

1  ((DZMDZ2)/0XIM((KXMKX21/2M.3»((0Z2*DZI)/0X2I) 

H(3,7)  =  (KY2*HM,7)»((DI2»DX3)/DZ2)*KYM(  (0X2*0X31/ 

1  DZ1MH(2,7M((KXMKX2)/2M.5*U01M0Z2I/DX3I»H(3,8M 
1  l(KX2*KXn/2)».5*UD2MDZ2>/DI2MH(3,6))/UKY2/DZ2M 
1  (DX2+DX3MKYM  1(0X2*0X31  /D2l)*(  (KX1+KX2I/210* 

1  ((DZ1*DZ2)/DX3)*((KX1*KX2I/2)*.5*<(0Z2*DZI>/DX2I) 

C  RON  TNO  -  INTERIOR  NODES 

DO  1700  C=2 , 5 

1700  H(2,CM  l(KXH(H(2,C*l!*H(2,C-l))  )/Hj*KY1*(H(3,C)*H(1 ,0) » 
l  H5) / ( (2*KX1)  /M5*2*KYI*(15) 

DO  1800  C=B, 1 1 

1800  H(2,CM((YXM(H(2,C-l)*H<2,C*ll))/H6*KYl»(H(3,C)*H(l,C)M 

1  t16)/((2*KXl)/l1S*2*KYl*(16) 

H(2,6)-<YXM(DZ1/DX2MHI2,7,+KXM!0ZMDXI)»H(2,51^YM 
1 < 1DX1+DX2) /  12*DZl) ) *<H(1 ,6) +H(3,6)  >)/<KXl*(DZ 1/0X2* 

Figure  A5.  (Sheet  6  of  8) 


A 1  3 


1  OZ 1/0X1) *KY1 *1(0X1 *0X2) /DZ1 ) ) 


H(2,7)=(KXI»(DZ1/0X3)*H(2,8)*KX1*(DZI/0X2)»H(2,6)*KY1» 

1  (IQX2*DX3)/(2*DZI))»(H(1,7)*H(3,7>))/(KXI*(DZ1/0X3* 

1  0Z1/0X2)*KY1*(  (0X2*0X31 /DID) 

C  SEARCH  FOR  HA XI HUH  RESIDUAL  HR 

HR=0.0 

DO  2200  R=l,8 
00  2200  C=1,12 
HRN=ABS (H (R , 0 -OLOH (R , C ) ) 

IF(flfiN.GT.HR)  HR=HRN 
2200  CONTINUE 

C  RETURN  FOR  NEXT  ITERATION 
IF(NR.6T.TQU  60  TO  300 

C  PRINT  COHPLETION  HESSAGE 

KRITE (20,950)  ITT,  NR 
WRITE (* ,950)  ITT, NR 

m  FORHATf/'  66UUT ION  CONPLITI  AFTER  ',15,'  ITERATIONS' 

1/'  NAXIHUH  RESIDUAL  =  79.5) 

60  TO  2450 

C  TERN I NAT I ON  FOR  NAXIHUH  TRIES 

2400  HRITE(20,952)  ITT,  NR 
WRITE (» ,952)  ITT, HR 

952  FORHAT (/ ’  NO  SOLUTION  AFTER  NAXIHUH  OF  ',15,*  ITERATIONS' 
1  /'  NAXIHUH  RESIDUAL  =  79.5) 

60  TO  2450 

C  CALCULATE  FINAL  HEADS  AND  GRADIENTS 

2450  RHO 1= (H(3,7)-H(l,7)) 

RHD2=(H(5,7)-H(3,7I) 

GRI-RHD1/Z1 

6R2=RHD2/Z2 

C  NRITE  FINAL  HEADS  TO  FILE 

HRITE(20,940I  UH(R,C),C=l,12),R=l,fl) 

940  FORMAT ( / '  FINAL  HEADS', /,12(1X,FL.;)) 


C  PRINT  RESIDUAL  HEAOS  AND  GRADIENTS 

WRITE (20,960)  RHD1,6R1,RHD2,6R2 
«RITE(«,960)  RH0l,GRl,RHD2,6R2 

960  FORHAK/'  RESIDUAL  HEAD  AT  RASE  OF  TOP  BLANKET  :  76. 2 
I/'  GRADIENT  THROUGH  TOP  BLANKET  =  76.2 
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21/'  HEAD  LOSS  THROUGH  MIDDLE  STRATUM  =  76.2 

If  GRADIENT  THROUGH  MIDDLE  STRATUM  =  76.2) 

C  CLOSE  FILES  AND  CHECK  FOR  MEN  PROBLEM 

REMIND  10 
REMIND  20 
CLOSE  (10) 

CLOSE  (20) 

MRITE  (*,970)  FILEOUT 
970  FORMAT!/'  OUTPUT  SAVED  IN  FILE  *,A) 

MRITE (* ,980) 

960  FORMAT!/'  MANT  TO  RUN  ANOTHER  PROBLEM  ?  Y  OR  N) 
READ!*,’ (A) ’I  ANS 

IF( (ANS.Efi,  Tl.OR. (ANS.EQ, 'y'D  GOTO  100 

STOP 

END 

C  INTRODUCTION  SUBROUTINE 

SUBROUTINE  INTRO 
MRITE (• ,20) 

20  FORMAT!/, 

I'  PROGRAM  LEVEE JL  —  7 

2'  UNOERSEEPAGE  ANALYSIS  FOR  3-LAYER  PROFILES'// 

3'  Mrltttn  by  Thaeas  F.  Naif  I  and  H.  A.  Al-Moussa*ti  7 
4'  Michigan  State  University 7/ 

5'  For  the  U.  S.  Arey  Corps  of  Engineers'/ 

6'  Release  1.0  September  1987') 

MRITE(*,30) 

30  FORMAT (///, 


1/'  -- 

/ 

2/'  . 

mm 

mx  xxm  tm 

3/'  . 

x 

XI  XXX  XX 

4/'  . 

mm 

mmxiixxxxxx 

5/'  . 

i 

XXI  1X1  1X1  IX 

6/'  . 

7/'  -- 
8//) 

mm 

xxxxmi  ixxxxxx 

RETURN 

END 
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APPENDIX  B:  COMPUTER  PROGRAM  LEVEE IRR:  UNDERSEEPAGE 
ANALYSIS  FOR  FOUNDATIONS  CHARACTERIZED  BY  TWO 
LAYERS  OF  IRREGULAR  SHAPE 


1.  The  program  LEVEEIRR  performs  underseepage  analysis  for  levee  foun¬ 
dations  consisting  of  two  layers  of  variable  thickness  with  irregular  bound¬ 
aries.  The  program  was  written  in  FORTRAN77  and  runs  on  IBM  PCs  or  compatible 
microcomputers  using  the  MS-DOS  operating  system. 

2.  Input  to  the  program  is  from  a  data  file  without  line  numbers.  The 
geometry  of  the  two  foundation  layers  is  described  by  dividing  the  foundation 
into  a  series  of  horizontal  segments.  At  x-coo-dinates  of  the  lines  defining 
the  segments,  the  elevations  of  the  top  of  the  top  blanket,  top  of  the  sub¬ 
stratum,  and  bottom  of  the  substratum  are  specified.  Up  to  30  segments  can  be 
specified,  and  it  is  recommended  that  minimum  of  3  segments  be  specified  on 
each  side  of  the  levee.  Node  points  at  which  the  flow  equation  is  solved  are 
generated  at  the  specified  top-of-substratum  coordinates  at  a  nine  additional 
points  between  each  set  of  adjacent  specified  points.  The  input  file  can  be 
created  using  any  word  processing  or  text-editing  program  that  produces  a 
standard  ASCII  file.  An  example  input  file  is  shown  in  Figure  Bl.  The  cross 
section  represented  by  the  example  is  shown  in  Figure  B2. 

3.  The  program  is  executed  by  typing  the  command  LEVEEIRR  with  the  file 
LEVEEIRR.EXE  resident  on  the  default  drive.  The  program  will  then  ask  for  the 
name  of  the  Input  file  and  the  name  of  the  output  file.  If  the  output  file 
already  exists,  it  will  be  written  over;  otherwise,  it  will  be  created.  A 
sample  run  of  LEVEEIRR  is  shown  in  Figure  B3.  An  example  output  file  is  shown 
in  Figure  BA.  A  program  source  liscing  of  LEVEEIRR  is  shown  in  Figure  B5 . 

A.  The  program  provides  output  to  two  devices:  the  console  (screen) 
and  the  output  file.  Output  to  the  screen  includes  the  number  of  iterations, 
the  residual  head  at  the  levee  toe,  and  the  gradient  at  the  levee  toe.  If  the 
user  desires  more  detailed  output,  the  output  file  can  be  displayed  or  printed 
using  any  word  processing  or  text  editing  program  that  works  with  standard 
ASCII  files.  The  output  file  contains  the  values  of  the  x-coordinates , 
y-coordinates ,  layer  thicknesses,  heads,  and  gradients  at  each  node  point. 


Bl 


EXAMPLE  INPUT  FILE  FOR  LEVEEIRR 


The  section  analyzed  using  this  data  file  is  illustrated  in  Figure  B2. 
Although  a  minimum  of  three  segments  on  each  side  of  the  levee  are  recom¬ 
mended,  only  four  segments  total  are  used  in  this  example  for  simplicity. 


Values  Variable  Names 


IRREGULAR  FOUNDATION  TEST  PROBLEM 

TITLE 

175.0  160.0 

HR  HL 

4  2 

NSEGS,  LSEG 

.2  .0002 

KF,  KB 

1000  .003 

TRIES,  T0L 

750.0  60.0 

140.0 

158.0 

X(l),  Y1 (1) ,  Y2 ( 1 ) 

,  Y3(l) 

1750.0  60.0 

140.0 

160.0 

X(2) ,  Y1 (2) ,  Y2 (2) 

,  Y3(2) 

1900.0  60.0 

140.0 

160.0 

X (3) ,  Y1 (3) ,  Y2(3) 

,  Y3 (3) 

2400.0  60.0 

120.0 

158.33 

X(4),  Y1 (4) ,  Y2(4) 

,  Y3(4) 

4900.0  60.0 

140.0 

150.0 

X(LAST) ,  Y 1 (LAST) 

Y2(LAST) 

Y3 (LAST) 


The  variables  are  defined  as  follows: 

TITLE  is  descriptive  title  for  the  problem  (80  characters  maximum). 

HR  is  the  elevation  of  water  on  the  riverside  of  the  levee.  Actual  elevations 
are  usually  most  convenient. 

HL  is  the  elevation  of  water  on  the  landside  of  the  levee  or  the  elevation  of 
landside  ground  if  no  water  is  present. 

NSEGS  is  the  number  of  segments  used  to  describe  the  problem  geometry. 

LSEG  is  the  segment  number  under  the  levee.  One  and  only  one  segment  must  be 
under  the  levee.  For  the  example  shown,  LSEG  equals  2. 

Figure  Bl.  Example  input  file,  program  LEVEEIRR  (Continued) 
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KF  is  the  horizontal  permeability  of  the  substratum,  in  feet  per  minute  or 
other  consistent  units. 

KB  is  the  vertical  permeability  of  the  top  blanket,  in  feet  per  minute  or 
other  consistent  units. 

TRIES  is  the  maximum  number  of  iterations  that  the  program  will  be  allowed  to 
make.  It  is  provided  to  stop  the  program  in  the  event  of  non-converging  solu¬ 
tions.  Generally  the  program  should  terminate  because  TOL  is  reached  before 
TRIES.  A  value  of  500  is  suggested. 

TOL  is  the  maximum  tolerance,  or  maximum  difference  in  head  at  any  node 
between  successive  iterations.  The  iteration  process  will  stop  when  the  maxi¬ 
mum  residual  is  smaller  than  TOL.  A  value  of  0.001  is  suggested. 

X(l),  Yl(l),  Y2(l),  Y3(l)  are  the  x-coordinate ,  the  y-coordinate  of  the  base 
of  the  substratum,  the  y-coordinate  of  the  base  of  the  top  blanket,  and  the 
y-coordinate  of  the  ground  surface,  respectively,  at  the  riverside  of  the 
first  segment.  This  line  is  repeated  for  the  riverside  of  each  segment. 

X(LAST) ,  Y 1 (LAST) ,  Y2(LAST),  Y3(LAST)  are  the  x-coordinate,  the  y-coordinate 
of  the  base  of  the  substratum,  the  y-coordinate  of  the  base  of  the  top 
blanket,  and  the  y-coordinate  of  the  ground  surface,  respectively,  at  the  land 
side  of  the  last  segment. 


Figure  Bl.  (Concluded) 
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Riverside  Landside 
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Figure  B2.  Definition  of  variables,  LEVEEIRR 


c> 

C  >LEVEE I RR 


PROGRAM  UEVEE I RR  - 

UNOERSEEPAOE  ANALYSIS  FOR  IRREGULAR  PROFILES 

Written  by  Thomas  F.  Wolf  I  and  H.  A.  Al-Mauasawi 
Michigan  State  University 

For  the  U.  S*  Army  Carps  o/  Engineers  Release  1.0  September  1907 


XXXXXX 

X 

XXXXXX 

X 

XXXXXX 


xxxx  XXXXX  XXXX 
XX  XXX  XX 
xxxxxxxi 1XXXXXX 
xxr  ixx  ixz  ix 

XXXXXXXI  IXXXXXX 


ENTER  INPUT  FILE  NAME 
DATAIRR 

ENTER  OUTPUT  FILE  NAME 
OUTIRR 

IRREGULAR  FOUNDATION  TEST  PROBLEM 


INPUT  OATA 

RIVERSIDE  HEAD 
173.00 

NO  OF  SEGMENTS 

4 

SUBSTRATUM  KF 

.200000 

MAX  ITERATIONS 
lOOO 


LAND8IDE  HEAD 
160.00 

SEGMENT  AT  LEVEE 
2 

BLANKET  KB 

.000200 

TOLERANCE 

.003000 


X  DISTANCE 
750.00 
1750.00 
1900.00 
2400.00 
4900.00 


BASE  OF  SUBSTRATUM 
60.00 
60.00 
60.00 
60.00 
60.00 


BASE  OF  BLANKET 
140.00 
140.00 
140.00 
120,00 
140.00 


TOP  OF  BLANKET 
150.00 
160.00 
160.00 
150.33 
150.00 


JUST  A  MOMENT,  I  AM  THINKING 


THE  SOLUTION  IS  COMPLETE 
TOTAL  ITTERATIONS*  298 
MAXIMUM  RESIDUAL-  .0030 


HEAD  AT  LEVEE  TOE  -  7. 38 

MAXIMUM  HEAD  -  7.30  AT  X  -  1900.00 

GRADIENT  AT  LEVEE  TOE  -  .37 

MAXIMUM  GRAOIENT  •  .37  AT  X  -  1900.00 

DETAILED  OUTPUT  SAVED  IN  FILE  OUTIRR 


WANT  TO  RUN  ANOTHER  PRQ8LEM  7  Y  OR  N 
N 

Stop  -  Program  terminated. 

C> 


Figure  R3.  Example  run  for  program  LEVEE1RR 
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EXAMPLE  OUTPUT  FILE  FOR  LEVEEIRR 


IRREGULAR  FOUNDATION  TEST  PROBLEM 

INPUT  DATA 

RIVERSIDE  HEAD  LANDSIDE  HEAD 

175.00  160.00 

NO  OF  SEGMENTS  SEGMENT  AT  LEVEE 

A  2 

SUBSTRATUM  KF  BLANKET  KB 

.200000  .000200 

MAX  ITERATIONS  TOLERANCE 

1000  .003000 

X  DISTANCE  BASE  OF  SUBSTRATUM  BASE  OF  BLANKET  TOP  OF  BLANKET 


750.00 

60.00 

1*0.00 

158.00 

1750.00 

60.00 

1*0.00 

160.00 

1900.00 

60.00 

1*0.00 

160.00 

2*00.00 

60.00 

120.00 

158.33 

*900.00 

60.00 

1*0.00 

150.00 

X  VALUES 

750.00  850.00  950.00  1050.00  1150.00  1250.00  1350.00  1*50.00  1550.00  1650.00 
1750.00  1765.00  1780.00  1795.00  1810.00  1825.00  18*0.00  1855.00  1870.00  1885.00 
1900.00  1950.00  2000.00  2050.00  2100.00  2150.00  2200.00  2250.00  2300.00  2350.00 
2*00.00  2650.00  2900.00  3150.00  3*00.00  3650.00  3900.00  *150.00  **00.00  *650.00 
*900.00 


D  VALUES 


80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

80.00 

78.00 

76.00 

7*  .00 

72.00 

70.00 

68.00 

66.00 

6*  .00 

62.00 

60.00 

62.00 

6*. 00 

66.00 

68.00 

70.00 

72.00 

7*  .00 

76.00 

78.00 

80.00 

Z  VALUES 

18.00 

18.20 

18. *0 

18.60 

18.80 

19.00 

19.20 

19.  *0 

19.60 

19.80 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

21.83 

23.67 

25.50 

27.33 

29.17 

31,00 

32.83 

3*. 66 

36.50 

38.33 

35.50 

32.66 

29.83 

27.00 

2*  .  17 

21.33 

18.50 

15.67 

12.83 

10.00 

READS  AT  TOP  OF  BLANKET 

175.00  175.00  175.00  175.00  175.00  175.00  175.00  175.00  175.00  175.00 


Figure  B4.  Example  output  file,  program  LEVEEIRR  (Continued') 


B6 


175.00  173.50  172 
160.00  160.00  160 
160.00  160.00  160 
160.00 


00 

170 

50 

169 

00 

00 

160 

00 

160 

00 

00 

160 

00 

160 

00 

167.50  166.00  164 . 50  163.00  161.50 
160.00  160.00  160.00  160.00  160.00 
160.00  160.00  160.00  160.00  160.00 


INITIAL  HEADS 


175 

.00 

174 

.50 

174 

.00 

173. 

.50 

173 

.00 

170. 

.00 

169 

.50 

169 

.00 

168. 

.50 

168 

.00 

165. 

.00 

164 

.50 

164. 

.00 

163. 

,50 

163 

.00 

160 

.00 

159 

.50 

159 

.00 

158 

.50 

158 

.00 

160.00 

THE  SOLUTION  IS  COMPLETE 
TOTAL  ITTERATIONS-  298 
MAXIMUM  RESIDUAL-  .0030 

X  DISTANCE 
PIEZ  EL 
RESIDUAL  HEAD 
GRADIENT 


172.50  172.00  171.50  171.00  170.50 

167.50  167.00  166.50  166.00  165.50 

162.50  162.00  161.50  161.00  160.50 

157.50  157.00  156.50  156.00  155.50 


750.0 

850.0 

950.0 

1050.0 

1150.0 

175.00 

174.40 

173.80 

173.19 

172.57 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

1750.0 

1765.0 

1780.0 

1795.0 

1810.0 

168.41 

168.29 

168.18 

168.07 

167.96 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

1900.0 

1950.0 

2000.0 

2050.0 

2100.0 

167.38 

167.09 

166.82 

166.56 

166.32 

7.38 

7.09 

6.82 

6,56 

6.32 

.37 

.32 

.29 

.26 

.23 

2400.0 

2650.0 

2900.0 

3150.0 

3400.0 

165.04 

164.16 

163.41 

162.76 

162.19 

5.04 

4.16 

3.41 

2.76 

2.19 

.  13 

.12 

.10 

.09 

.08 

HEAD  AT 

LEVEE  TOE  - 

7.38 

MAXIMUM 

HEAD  - 

7.38  AT 

X  -  19! 

1250.0 

1350.0 

1450.0 

1550.0 

1650.0 

171.94 

171.28 

170.61 

169.91 

169.17 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

182S.0 

1840.0 

1855.0 

1870.0 

1885.0 

167.85 

167.75 

167.65 

167.56 

167.46 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

2150.0 

2200.0 

2250.0 

2300.0 

2350.0 

166.08 

165.85 

165.63 

165.43 

165.23 

6.08 

5.85 

5.63 

5.43 

5.23 

.21 

.  19 

.  17 

.16 

.  14 

3650.0 

3900.0 

4150.0 

4400.0 

4650.0 

161.70 

161.28 

160.91 

160.58 

160.28 

1.70 

1.28 

,  91 

.58 

.28 

.07 

.06 

.05 

.04 

.02 

.00 


GRADIENT  AT  LEVEE  TOE  -  .37 

MAXIMUM  GRADIENT  -  .37  AT  X  '  1900.00 


Figure  B4.  (Concluded) 
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PROGRAM  IEVEE1RR 


ZMMMtmttttnnutiMinmmmmmmtMMfmmitiu 

c 

C  THIS  PROGRAM  HAS  WRITTEN  BY: 

C 

C  DR.  THOMAS  F.  WOLFF;  ASSISTANT  PROFESSOR 

C  AND 

C  HASSAN  M.  ALNOUSSANI;  PH.D  CANDIDATE 

C 

C  US INS  THE  A.  H.  CASE  CENTER 

C  COMPUTER  AIDED  ENSINEE1NS  FACILITY  AT 

C  MICHIGAN  STATE  UNIVERSITY  ,  1987 

C 

zmtntmtmmmtmmtmmmummummmmmmm 

c 

C  THIS  PROGRAM  USES  THE  FINITE  DIFFERENCE  METHOD  TO  ANALYZE 
C  SEEPAGE  THROUGH  TWO  SOIL  LAYERS  BENEATH  A  LEVEE. 

C  THE  TOP  SOIL  LAYER  HAS  PERMEABILITY  KB  AND  THICKNESS  Z 
C  THE  BOTTOM  LAYER  HAS  PERMEABILITY  KF  AND  THICKNESS  D 
C  Z  AND  D  ARE  FUNCTIONS  OF  THE  DISTANCE  It  FROM  THE  LEFT. 

C 

C 


c 

r 

VARIABLES: 

l 

c 

A 

Segaent  Nuaber  under  levee 

c 

CLEFT  (0 

Left  side  head  coefficient 

c 

CRIGHT(C) 

Right  side  head  coefficient 

c 

CLEV1 

Coluan  no.  at  riverside  levee  toe 

c 

CLEV2 

Coluan  no.  at  landside  levee  toe 

c 

CTOP(c) 

Top  side  head  coefficient 

c 

DX(C) 

Differential  foundation  width  at  node  C 

c 

HL 

Landside  water  elevation 

c 

HR 

Riverside  water  elevation 

c 

H(R,C) 

Head  (row  1  is  ground,  row  2  is  base  of  blanket) 

c 

ITT 

Iteration 

c 

KB 

Peraeability  of  blanket 

c 

KF 

Peraeability  of  substratua 

c 

LENGTH 

Segaent  Length 

c 

MR 

Maxiaua  residual 

c 

N 

Nuaber  of  foundation  segaents 

c 

NN 

Node  subscript 

c 

OLDHIR.C) 

Head  fro*  previous  iteration 

c 

TOL 

Tolerance 

c 

TRIES 

Maxiaua  nuaber  of  iterations 

c 

HR,C> 

X  coordinate  at  coluan  C 

c 

YKC) 

Y  coordinate  at  base  of  substratua 

c 

Y2(C) 

Y  coordinate  at  base  of  blanket 

c 

Y3(C) 

Y  coordinate  at  ground  surface 

c 

Z 

Top  blanket  thickness 

Figure  B5.  Listing  of  program  LEVEETRR  (Sheet  1  of  8) 
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C  DECLARE  VARIABLES 


CHARACTER  ANS#1,  FILEIN«Z>4,  FILE0UT*64,  TITLE»80 
INTEGER  A,C, Cl , C2 , CLEV 1 , CLEV2 , I TT , J ,N ,HN ,NN 10 . TRIES 
REAL  HR, HL , H (2 , 300) ,RHHAX,DH,DN1,0LDH(300) 

REAL  X (300) ,DX(300) ,DXCTR(300) ,Y1 (300) ,Y2(300) ,Y3(300> 
REAL  Z ( 300) , D (300) ,RESHD (300) ,8RAD (300) , SRNAJ 
REAL  CLEFT (300) , CR1SHT 1300) ,CT0P (300) ,CSUH (300) 

REAL  KB,KF,«RN,NR,T0L 
REAL  FC,FJ,FNN, LENGTH 
REAL  XHMAX ,  X6RHAX 

C  OISPLAY  INTRODUCTION 
CALL  INTRO 

C  OPEN  DATA  FILES 

100  WRITE  (*,*>  'ENTER  INPUT  FILE  NAME' 

READ  <V<A)')  FILEIN 

OPEN  (UNIT-10,  FILE=FILEIN,  ERR=100,  STATUS*  OLD') 
REMIND  (10) 

WRITE  (♦,«)  ENTER  OUTPUT  FILE  NANE' 

READ ( * , ' (A) ' )  FILEOUT 
OPEN  (20,FILE*FILE0UT) 

REMIND  (20) 

C  READ  DATA 

READ ( 10 , ' (A) ' ) TITLE 

READ ( 10, *)  HR,HL,N,A,KF,KB,TRIES,TOL 

NN=N*10 

NN1C=NN+10 

READ ( 10, *)  (X(J) ,Y1 (J) ,Y2(J> ,Y3(J) ,J=10,WW10,10> 

C  PRINT  INPUT  DATA  AND  MRITE  TO  FILE 

WRI TE ( 20 , * )  TITLE 
NRITEU,*) 

MRITE(*,«)  TITLE 

WRITE420 ,910) 

MR1TE(»,910) 

910  FORMAT!/ '  INPUT  DATA  ') 

WRI TE (20,932) 

MRITE(«,912) 

912  FORMAT)/,'  RIVERSIDE  HEAD  LANDSIDE  HEAD') 
NfiITE(20,914)  HR ,HL 
WRITE («,914>  HR , HL 
914  FORMAT  (1X,F  10.2,10)1,  FI  0.2) 

WRITE (20,916) 

NRITE<»,916) 

Figure  B5.  (Sheet  2  of  8) 
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SEGMENT  AT  LEVEE') 


916  FORMAT(/,‘  NQ  OF  SEGMENTS 
NRITE(2Q,918)  N.A 
WRITE!*, 918)  N,A 
918  FORMAT  <tX,I3,17X,I3) 

WRITE (20,920) 

WRITE!*, 920) 

920  FORMAT </,'  SUBSTRATUM  KF  BLANKET  KB') 

WRITE (20,922)  KF,KB 
WRITE (*,922)  KF, KB 
922  FORMAT I1X,F12.6,SX,F12.6> 

WRITE (20,924) 

WRITE!*, 924) 

924  FORMAT!/, '  MAX  ITERATIONS  TOLERANCE') 

WRITE (20,926)  TRIES, TOL 
WRITE!*, 926)  TRIES, TOL 
926  FORMAT  (1X,I5,15X,FJ0.6>) 

WRITE (20, 928) 

WRITE!*, 928) 

928  FORMAT!/,'  X  DISTANCE  BASE  OF  SUBSTRATUM  BASE  OF  BLANKET 
1T0P  OF  BLANKET') 

WRJTE<20,930)  !X(J> ,Y1(J) ,Y2!J) ,¥3U> ,M0,MN10,10) 
WRITE<»,930)  (X(J)  ,Y1U)  ,Y2(J!  ,Y3(J) , J=10,MN10, 10) 

930  F0RMAT(1X,4!F10.2,6X>) 


C  CALCULATE  THE  X  VALUES  FOR  ALL  NODES 

DO  200  CS10,NN,10 
LENSTH=X(C*10)-X(C> 

DX(C)*LENGTH/10.0 
DO  200  J=1 ,9 
DX(C*J)=DX(C) 

X(C*J)SX(C+J-1)*DX(C*J-1) 

200  CONTINUE 

C  WRITE  THE  X  VALUES  TO  FILE 

WRITE  (20,935)  <X(C),  C=10,NN*10) 

935  FORMAT  !/,'  X  VALUES' ,/, 10(1 XF7. 21 ) 

C  CALCULATE  THE  D,  Z,  AND  Y3  VALUES 
DO  350  C= 10 , NN , 10 
D(C)=Y2(C)-Y1 (C) 

Z(C)*Y3(C)-Y2(C) 

DO  340  J= 1 ,9 
FJ  =FLOAT (J) /10.0 

D(C*J)=Y2(C)-Y1(C)*FJ*(Y2(C*10)-Y2(C))-FJ*(Y1(C*!0)-Y1(C>) 

Z(C+J)=Y3(C)-Y2(C)*FJ*(Y3(CM0)-Y3(C)!-FJ*(Y2(CM0)-Y2(C)I 

Y3(C*J)=Y3(C)*FJ»(Y3IC*10)-Y3(CI) 

340  CONTINUE 
350  CONTINUE 
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D(NN10)=Y2(NN101-Y1(NN10) 

Z(NN10)=Y3(NN10)-Y2(NN10) 

C  WRITE  THE  D  AND  Z  VALUES  TO  FILE 

WR1TE<20,940! 

940  FORMAT (/, ’  D  VALUES' ) 

WRITE < 20, 941 2  1510,  C=10,NN+10) 

941  FORMAT (1041XF6.2) ) 

WRI TE (20 , 942) 

942  FORMAT ( / , '  Z  VALUES' ) 

WRI TE (20,941)  <Z(C),  C-iO.NN+lO) 

C  FIND  THE  COLUMNS  UNDER  THE  LEVEE  --  CLEV1  TO  CLEV2 

CLEV1  «  A*10 
CLEV2  =  CLEV1  ♦  10 

C  CALCUALATE  HEAD  CHANGES  AT  LEVEE/BLANKET  INTERFACE 

1F(YJ(CLEV2).GT,HL)  THEN 
DH= (HR-Y3 (CLEV2) 1/10.0 
ELSE 

DH= (HR-HL 1/10.0 
ENDIF 

C  CALCULATE  THE  HEAD  AT  TOP  OF  BLANKET 

DO  400  OlO.NNIQ 

C  RIVERSIDE 

IF(C.LE.CLEVl)  THEN 
H<1 ,C)=HR 

C  LEVEE 

ELSE  IF ( (C.6T.CLEV1), AND. (C.LT.CLEV2) )  THEN 
H(l,CI=H(l,C-l)-DH 

C  LANDS  IDE  -  WATER  ABOVE  GROUND 

ELSE1F I (C.6E.CLEV2) . AND. (HL.GT.YJ(C) ) I  THEN 
H(1,C)=HL 

C  LANDSIDE  -  GROUND  ABOVE  WATER 

ELSE 

Hd.Cl-YKC) 

ENDIF 

C  LANDSIDE  -  HIGH  GROUND  ABOVE  HR 

IF  I (C.GE.CLEV2) . AND. (HI1,C).6T.HR1)  THEN 
H(l.C)  =  Hll.C-D 
ENDIF 

40  CONTINUE 
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C  WRITE  HEADS  AT  TOP  OF  BLANKET  TO  FILE 
WRITE (20 ,?45) 

945  FORMAT!/, '  HEADS  AT  TOP  OF  BLANKET  ) 

WRI TE ( 20 , 94 1 )  (H<l,C>,  C=i0.NN10) 

C  INITIALIZE  HEADS  AT  BASE  OF  BLANKET 

H ( 2, 10) =HR 
H(2,NN10)=HL 
FNN=FLOAT (NN) 

DO  450  Csll,NN*9 
FC=FLOAT (C) 

H(2,C)SH(2,10M 'FC-10.0)/ (FNN-10.0) )*<HR-HL) 
450  CONTINUE 

C  WRITE  INITIAL  HEADS  AT  BASE  OF  BLANKET  TO  FILE 
WR1TE(20,950) 

950  FORMAT!/, '  INITIAL  HEADS') 

WR1TE(20,941)  (H!2,C),  C=10,NN10) 

i  Calculate  difference  coefficients 

DO  470  C=U,  NN+9 
CLEFT (C)=KF«D(C) /DX (C-l ) 

DXCTR!C)*(DX(C-1)+DX  101/2 
CTOP (C) =KB*DXCTR  <0 / Z (C) 

CRIGHT (C)=KF*D(C)/DX(C) 

CSUM  (C) =CLEFT  (0  +CRISHT  (0  +CTOP  (0 
470  CONTINUE 


[UtflHtHIHHHHtiltHiHHIHtHtid 

C  ITTtRATIQN  PROCESS 


C  INITIALIZE 
I  TT=0 

WRITE  <#,*i 

WRITE  (*,«)  'JUST  A  M0MES7,  I  AH  THINKING' 

C  BEGIN  NEXT  ITERATION 

500  ITT=1TTM 

C  WRITE!*, 9X0)  ITT 

C  960  FORMAT C 1 X , ' ITERATION  NO  ’,15) 

IF UTT  +  l.GT. TRIES)  GO  TO  2400 

C  RESET  OLD  HEADS 

DO  510  C=t 1 ,NN+9 
510  0LDH(C)=H12,C) 
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C  CALCULATE  HEADS  AT  RON  TNQ 


DO  520  C=Il,Nh+9 

H(2,C)  =  (H(2,C-1  >*CLEFT(C>-»-H(l  ,C:*iCT0PLCi^H(2,C+l)*CRIGHT  (0) 
1  /CSUMIC) 

520  CONTINUE 

C  SEARCH  FOR  MAXIMUM  RESIDUAL  -HR- 


MR=0.0 

DO  600  C=1 1 ,NN+9 
MRN=A8S  (H  (2,0  -DLDH  (C)  > 

IF (MRN.6T.MR)  MR=MRN 
600  CONTINUE 

C  RETURN  FOR  NEXT  ITERATION 
IF  (HR.  6T.  TOD  SO  TO  500 

C  NORMAL  TERMINATION 

WRITE (20,970) ITT, HR 
NRITE(*,970)  ITT, HR 

970  FORMAT (/ , '  THE  SOLUTION  IS  COMPLETE ' 
l/'  TOTAL  ITTERATIONS*  ',14 
2/'  MAXIMUM  RESIDUAL*  \F6.4) 

C  CALCULATE  FINAL  RESIDUAL  HEADS  AND  GRADIENTS 

620  DO  450,  C=CLEV2,NN*9 
RESHD(C)=H(2,C)-H(1,C) 

GRAD  (C)  =RESHD  (0 II  1C) 

650  CONTINUE 

C  CALCULATE  MAX  LANDS  IDE  RESIDUAL  HEAD 

RHHAX=RESHD(CLEV2) 

XHHAX=X (CLEV2) 

DO  660  C=CLEV2 , NN+9 
IF  (RESHD(C) . ST . RHMAX )  THEN 
RHHAX-RESHD(C) 

XHMAX=X (C ! 

ENOIF 
660  CONTINUE 

C  CALCULATE  HAX  LANDSIDE  6RADIENT 
SRMAX=6RAD (CLEV2) 

XSRHAX=X (CLEV2) 

DO  670  C=CLEV2,NN+9 
IF  (SRAD(C) .ST.GRHAX)  THEN 
6RHAX=6RAD(C) 

XGRHAX=X(C) 

ENDIF 
670  CONTINUE 
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C  WRITE  X,  PIEZ  EL,  HEAD,  AND  6RADIENT  TD  FILE 
WRITE  (20,975) 

975  FORMAT ( / , '  X  DISTANCE1,/,'  PIEZ  EL',/,'  RESIDUAL  HEAD1,/ 
1‘  GRADIENT',/) 

DO  700  1*1, N 
C1=I*10 
C2=I*10*9 

WRITE (20,980)  (HO,  C=C1,C2) 

WRITE(20,982)  (HI2,C),  C=C1,C2) 

NRITE(20,982t  (RESHD(C) ,  C=C1 ,C2) 

NRITE!20,982)  (GRAD(C),  C=C1,C2) 

700  CONTINUE 

980  FORMAT!/, 10(1XF7.1)) 

982  FORMAT (10(1XF7.2>) 

C  PRINT  TOE  AND  MAXIMUM  RESIDUAL  HEADS  AND  GRADIENTS 
WR1TE(20,990)  RESHD CCLEV2) ,  RHMAX,  XHMAX 
WRITE!*, 990)  RESHD (CLEV2) ,  RHMAX,  XHMAX 
990  FORMAT! /,'  HEAD  AT  LEVEE  TOE  *  ',F8.2,/ 

1'  MAXIMUM  HEAD  =  ',F8.2,'  AT  X  =  ‘  ,F8.2) 

WRITE (20,994)  GRAD1CLEV2) ,  SRMAX,  XBRMAX 
WRITE ( *,994)  6RAD (CLEV2) ,  SRMAX,  XSRMAX 

994  FORMAT!/,'  GRADIENT  AT  LEVEE  TOE  =  ' ,F8.2,/ 

1*  MAXIMUM  GRADIENT  *  ',F8.2,'  AT  X  *  ',F8.2) 

C  CLOSE  FILES  AND  CHECK  FOR  NEW  PROBLEM 
REWIND  10 
REWIND  20 
CLOSE  (10) 

CLOSE  (20) 

WRITE  (*,1000)  F1LE0UT 

1000  FORMAT!/'  DETAILED  OUTPUT  SAVED  IN  FILE  \A) 

WRITE!*, 1010) 

1010  FORMAT!/, •  WANT  TO  RUN  ANOTHER  PROBLEM  ?  V  OR  N  ') 
READ!*, ' (A) ‘)  ANS 
IF  (ANS.E8.T)  SOTO  100 
GO  TO  99 

C  TERMINATION  FOR  MAXIMUM  TRIES 

2400  1(RITE(20,995)  ITT, MR 
WRITE!*, 995)  ITT, MR 

995  FORMAT!//,  'SOLUTION  INCOMPLETE  STOPPED  BY  MAXIMUM 

1  POSSIBLE  NUMBER  OF  TRIES*  1 , 14, /10X, 'MAXIMUM  RESIDUAL* 

1  ',F6.3) 

GO  TO  620 

99  STOP 

END 

Figure  B5.  (Sheet  7  of  8) 
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C  INTRODUCTION  SUBROUTINE 
SUBROUTINE  INTRO 
NRITE<*,20) 

20  FORMAT  (/, 

r  PROGRAM  LEVEEIRR  —  7 
V  UNDERSEEPAGE  ANALYSIS  FOR  IRREGULAR  PROFILES'// 

3'  Written  by  Tho»as  F.  Nolff  and  H,  A.  Al-Houssa*i 7 
Michigan  State  University'// 

5'  For  the  U,  S.  Ar«y  Corps  of  Engineers' 


A‘  Release  1.0 
NRITEU.30) 
FORMAT!///, 

1/'  . 

September  1987') 

2/'  .  mm 

XXXX  XXXXX  XXXX 

3/'  .  X 

XX  XXI  XX 

!/'  .  mm 

xmxxxmxxxxx 

5/'  .  X 

XXI  1X1  1X1  IX 

a r  .  mm 

v  - 

8//) 

RETURN 

END 

mmxi  imm 

Figure  B5.  (Sheet  8  of  8) 
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APPENDIX  C:  COMPUTER  PROGRAM  LEVEECOR:  UNDERSEEPAGE 
ANALYSIS  AT  ANGLES  OR  "CORNERS"  IN  LEVEE  ALIGNMENT 


1.  The  program  LEVEECOR  performs  underseepage  analysis  at  bends  or 
corners  in  levee  alignment.  The  program  was  written  in  FORTRAN77  and  runs  on 
IBM  PCs  or  compatible  microcomputers  using  the  MSD-DOS  operating  system. 

2.  Input  to  the  program  is  from  a  15  line  data  file  without  line 
numbers.  This  file  can  be  created  using  any  word  processing  or  text-editing 
program  that  produces  a  standard  ASCII  file.  An  example  input  file  is  shown 
in  Figure  Cl.  Program  variables  are  defined  in  Figure  C2. 

3.  The  program  is  executed  by  typing  the  command  LEVEECOR  with  the  file 
LEVEECOR.EXE  resident  on  the  default  drive.  The  program  will  then  ask  for  the 
name  of  the  input  file  and  the  name  of  the  output  file.  If  the  output  file 
already  exists,  it  will  be  written  over;  otherwise,  it  will  be  created.  A 
sample  run  of  LEVEECOP  is  shown  in  Figure  C3.  An  example  output  file  is  shown 
in  Figure  C4.  A  program  source  listing  of  LEVEECOR  is  shown  in  Figure  C5. 

4.  The  program  provides  output  to  two  devices:  the  console  (screen) 
and  the  output  file.  Output  to  the  screen  includes  the  number  of  iterations, 
the  residual  head  at  the  levee  toe,  and  the  gradient  at  the  levee  toe.  If  the 
user  desires  more  detailed  output,  the  output  file  can  be  displayed  or  printed 
using  any  word  processing  or  text-editing  program  that  works  with  standard 
ASCII  files.  The  output  file  contains  the  approximate  distances  to  adjacent 
nodes,  equation  coefficients,  and  initial  and  final  heads  at  each  node  point. 
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EXAMPLE  INPUT  FILE  FOR  LEVEECOR 


Values 

TEST  PROBLEM 
12.58 
0 

5000 

650 

300 

80 

.089 

.0000908 

.0000908 

1100.0 

75 

.005 

2000 

14 


Variable  Names 
TITLE 
HR 
HL 
L1A 
LIB 
L2 

THETA 

KF 

KBL 

KBR 

RADIUS 

D 

TOL 

TRIES 

Z 


The  variables  are  defined  as  follows: 

TITLE  is  descriptive  title  for  the  problem  (80  characters  maximum) . 

HR  is  the  elevation  of  water  on  the  riverside  of  the  levee.  If  the  landside 
water  or  landside  ground  is  taken  as  the  datum,  HR  is  the  net  head  on  the 
levee . 

HL  is  the  elevation  of  water  on  the  landside  of  the  levee  or  the  elevation  of 
landside  ground  if  no  water  is  present.  This  level  is  often  taken  as  the 
datum,  in  which  case  HL  is  0.0. 

L1A  is  the  distance  from  the  riverside  levee  toe  to  the  river  or  source  of 
seepage  on  the  left  side  of  the  bend  if  one  is  facing  the  river. 

Figure  Cl.  Example  input  file,  program  LEVEECOR  (Continued) 
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LIB  is  the  distance  from  the  riverside  levee  to  to  the  river  or  source  of 
seepage  on  the  right  side  of  the  bend  if  one  is  facing  the  river. 


L2  is  the  base  width  of  the  levee. 

THETA  is  the  deflection  angle  between  two  levee  tangents  (analogous  to  the 
delta  angle  in  surveying) . 

KF  is  the  horizontal  permeability  of  the  substratum,  in  feet  per  minute  or 
other  consistent  units. 

KBL  is  the  upward  vertical  permeability  of  the  landside  top  blanket,  in  feet 
per  minute  or  other  consistent  units. 

KBR  is  the  downward  vertical  permeability  of  the  riverside  top  blanket,  in 
feet  per  minute  or  other  consistent  units. 

RADIUS  is  the  radius  of  curvature  of  the  landside  levee  toe,  in  feet  or  other 
consistent  units. 

D  is  the  thickness  of  the  substratum  in  feet  or  any  consistent  units. 

TRIES  is  the  maximum  number  of  iterations  that  the  program  will  be  allowed  to 
make.  It  is  provided  to  stop  the  program  in  the  event  of  non-converging 
solutions.  Generally  the  program  should  terminate  because  TOL  is  reached 
before  TRIES.  Because  LEVEECOR  may  converge  slowly,  a  value  of  1,000  is 
suggest -d . 

TOL  is  the  maximum  tolerance,  or  maximum  difference  in  head  at  any  node 
between  successive  iterations.  The  iteration  process  will  stop  when  the 
maximum  residual  is  smaller  than  TOL.  Because  LEVEECOR  may  converge  slowly,  a 
value  of  0.001  is  suggested. 

Z  is  the  total  thickness  of  the  top  stratum. 

Figure  Cl.  (Concluded) 
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c> 

OLEVEECOR 
PROGRAM  LEVEECOR 

SEEPAGE  ANALYSIS  AT  LEVEE  CORNERS 

ENTER  INPUT  FILE  NAME 
DATACOR 

ENTER  OUTPUT  FILE  NAME 
OUTCOR 


INPUT  PARAMETERS  OF  THE  PROGRAM 


TEST  PROBLEM 


HR  VALUE:  12.580 
HL  VALUE:  .000 

L1A  VALUE:  5000.000 
LIB  VALUE:  650.000 
L2  VALUE:  300.000 
THETA  VALUE:  80.000 
KF  VALUE: 

KBL  VALUE:  .000071 
KBR  VALUE:  .000091 
RADIuS  VALUE:  1100.000 
D  VALUE:  75.000 
TOL  VALUE:  .005000 
TOTAL  NUMBER  OF  TRIES:  2000 

thickness  value:  u.ooooqo 


ITERATION 

1 

HEAD  * 

6.05448 

MAI 

RES  = 

-5.54076  AT 

2, 

14 

ITERATION 

2 

HEAD  = 

5.95314 

MAI 

RES  = 

1.56634  AT 

5, 

9 

ITERATION 

3 

HEAD  * 

5.89575 

MAI 

RES  = 

.94499  AT 

5, 

7 

ITERATION 

A 

HEAD  = 

5.91939 

MAI 

RES  = 

.64477  AT 

5, 

6 

ITERATION 

5 

HEAD  = 

6.00586 

MAI 

RES  = 

.39485  AT 

5, 

6 

Figure  C3.  Example  run  for  program  LEVEECOR  (Continued) 
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ITERATION 

13 

HEAD  = 

6.34391 

HA!  RES  = 

-.04920  AT 

2,  10 

ITEKATION 

14 

HEAD  = 

6.3587! 

HA!  RES  = 

-.03832  AT 

2,  10 

ITERATION 

15 

HEAD  '= 

6.37394 

HA!  PES  = 

-.02976  AT 

2,  6 

ITERATION 

16 

HEAD  = 

6.38382 

HA!  RES  = 

-.02307  AT 

2,  6 

ITERATION 

17 

HEAD  = 

6.39376 

HA!  RES  * 

-.01784  AT 

2,  6 

ITERATION 

18 

HEAD  = 

6.40037 

HA!  RES  = 

-.01376  AT 

2,  10 

ITERATION 

19 

HEAD  * 

6.40693 

HA!  RES  * 

-.01058  AT 

2,  10 

ITERATION 

20 

HEAD  = 

6.41137 

HA!  RES  ■= 

-.00811  AT 

2,  10 

ITERATION 

21 

HEAD  ■ 

6.41577 

HA!  RES  ' 

-.00620  AT 

2,  6 

TOTAL  NUHBER  OF 

ITERATIONS 

21 

TEST  PROBLEM 


VALUE  OF  THE  6RA0IENT  AT  H<3,8)  *  .4583 

VALUE  OF  THE  H13,8)  «  6.4158 

Stop  -  Proqrii  ter«inated. 


C> 


Figure  C3.  (Concluded) 
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EXAMPLE  OUTPUT  fILE  FOR  LEVEECOR 


INPUT  PARAMETERS  OF  THE  PROGRAM 


TEST  PROBLEM 

HR  VALUE)  12.580 
HL  VALUE!  .000 
L1A  VALUE!  5000.000 
LIB  VALUE!  650.000 
L2  VALUE!  300.000 
THETA  VALUE!  80.000 
KF  VALUE!  .089000 
KBL  VALUE!  .000091 
KBR  VALUE!  .000091 
RADIUS  VALUE!  1100.000 
D  VALUE!  75.000 
TQL  VALUE!  .003000 
TOTAL  NUMBER  OF  TRIES!  2000 
THICKNESS  VALUE!  14.000000 
INITIAL  MATRIX 


.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6,290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6.290 

6,290 

6.290 

6.290 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

12.58 

DXLEFT  MATRIX 
.0000  1100. 

1100. 

1100. 

1100. 

1037. 

20.00 

20.00 

20.00 

20.00 

1037. 

1100. 

1100. 

1100. 

1100. 

.0000 

1100. 

1100. 

1100. 

1100. 

1100. 

20.00 

20.00 

20,00 

20.00 

1100. 

1100. 

1100. 

1100. 

1100. 

.0000 

1100. 

1100. 

1100. 

1100. 

1100. 

384.0 

384.0 

384.0 

384.0 

1100. 

1100. 

1100. 

1100. 

1100. 

.0000 

1100. 

1100. 

1100. 

1100. 

1100. 

488.7 

488.7 

488.7 

488.7 

1100. 

1100. 

1100. 

1100. 

1100. 

.0000 

1100. 

1100. 

1100. 

1100. 

1100. 

1172. 

981.7 

791.9 

602.1 

1100. 

1100. 

1100. 

1100. 

1100. 

.0000 

1100. 

1100. 

1100. 

1100. 

1100. 

1854. 

1475. 

1095. 

715.6 

1100. 

1100. 

1100. 

1100. 

1100. 

DXRIGHT  MATRIX 

Figure  C4.  Example  output  file,  program  LEVEECOR 
(Sheet  1  of  5) 
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20.00 


1100. 

20.00 

1100. 

20.00 

1100. 

J84.0 

1100. 

488.7 

1100. 

1172. 

1100. 

1854. 


1100. 

1037. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 


1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 


11C0. 

.too. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 


1037.  20.00 

1100.  1100. 

1100.  20.00 

1100.  1100. 

1100.  384.0 

1100.  1100. 

1100.  488.7 

1100.  1100. 

1100.  602.1 

1100.  1100. 

1100.  715.6 

1100.  1100. 


20.00 
.0000 
20.00  20.00 

.0000 

384.0  384.0 

.0000 

488.7  488.7 

.0000 

791.9  981.7 

.0000 

1095.  1475. 

.0000 


DYTOP  MATRIX 


75.00 

75.00 

75.00 

75.00 

75.00 

97.91 

97.91 

97.91 

97.91 

75.00 

75.00 

75.00 

75.00 

75.00 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

2500. 

2500. 

2500. 

2500. 

2500. 

2500. 

4456. 

3369. 

325.0 

325.0 

325.0 

325.0 

325.0 

325.0 

2500. 

2S00. 

2500. 

2500. 

2500. 

2500. 

4456. 

3369. 

325.0 

325.0 

325.0 

325.0 

325.0 

325.0 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

97.91 

1100. 

300.0 

3913. 

3913. 

.0000 


DYBOT  MATRIX 

.0000  .0000  .0000 

.0000  .0000  .0000 

75.00  75.00  75.00 

97.91  97.91  75.00 

1100.  1100.  1100. 

1100.  1100.  1(00. 

300.0  300.0  300.0 

300.0  300.0  300.0 

2500.  2500.  2500. 

3369.  325.0  325.0 

2500.  2500.  2500. 

3369.  325.0  325.0 


.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

75.00 

75.00 

97.91 

97.91 

97.91 

75.00 

75.00 

75.00 

75.00 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

1100. 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

2500. 

2500. 

2500. 

4456. 

3913. 

325.0 

325,0 

325.0 

325.0 

2500. 

2500. 

2500. 

4456. 

3913. 

325.0 

325.0 

325.0 

325.0 

DYCTR  MATRIX 
37.50  37.50 

48.95  48.95 

587.5  587.5 

599.0  599.0 

550.0  550.0 

550.0  550.0 

1250.  1250. 

1684.  162.5 

2500.  2500. 

3369.  325.0 


37.50  37.50  37.50 

37.50  37.50  37.50 

587.5  587.5  587.5 

587.5  587.5  587.5 

550.0  550.0  550.0 

550.0  550.0  550.0 

1250.  1250.  1250. 

162.5  162.5  162.5 

2500.  2500.  2500. 

325.0  325.0  325.0 


48.95  48.95  48.95 

37.50  37.50 

599.0  599.0  599.0 

587.5  587.5 

550.0  550.0  550.0 

550.0  550.0 

1250.  2228.  1956. 

162.5  162.5 

2500.  4456.  3913. 

325.0  325.0 


Figure  C4.  (Sheet  2  of  5) 
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1250. 

1250. 

1250. 

1250. 

1250. 

1250. 

2228. 

1956. 

1684. 

162.5 

162.5 

162.5 

162.5 

162.5 

162.5 

DXCTR  MATRIX 

550.0  1100. 

1100. 

1100. 

1069. 

528.5 

20.00 

20.00 

20.00 

528.5 

1069. 

1100. 

1100. 

1100. 

550.0 

550.0 

1100. 

1100. 

1100. 

1100. 

560.0 

20.00 

20.00 

20.00 

560.0 

1100. 

1100. 

1100. 

1100. 

550.0 

550.0 

1100. 

1100. 

1100. 

1100. 

742.0 

384.0 

384.0 

384.0 

742.0 

1100. 

1100. 

1100. 

1100. 

550.0 

550.0 

1100. 

1100. 

1100. 

1100. 

794.3 

488.7 

488.7 

488.7 

794.3 

1100. 

1100. 

1100. 

1100. 

550.0 

550.0 

1100. 

1100. 

1100. 

1100. 

851.1 

981,7 

981.7 

981.7 

851.1 

1100. 

1100. 

1100. 

1100. 

550.0 

550.0 

1100. 

1100. 

1100. 

1100. 

907.8 

1475. 

1475. 

1473. 

907.8 

1100. 

1100. 

1100. 

1100. 

550.0 

CLT  MATRIX 

.0000 

.3000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

3.365 

3.565 

3.565 

3.565 

3.635 

199.9 

199.9 

199.9 

199.9 

3.565 

3.565 

3.565 

3.565 

3.565 

.0000 

4.248 

4.248 

4.248 

4.248 

4.248 

12.17 

12.17 

12.17 

12.17 

4.248 

4.248 

4.248 

4.248 

4.248 

.0000 

8.493 

8.495 

8.495 

8.495 

8.495 

32.48 

28.77 

25.06 

4.268 

1.896 

1.896 

1.896 

1.896 

1.896 

.0000 

13.17 

13.17 

15.17 

15.17 

15.17 

25.39 

26.60 

28.39 

3.603 

1.972 

1.972 

1.972 

1.972 

1.972 

.0000 

.0000 

.0000 

.0000 

.0000 

.0006 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

CT  MATRIX 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

3.338 

6.675 

6.675 

6.675 

6.675 

3.398 

.1214 

.1214 

.1214 

3.398 

6.675 

6.675 

6.675 

6.675 

3.338 

12.24 

24.48 

24.48 

24.48 

24.48 

16.51 

8.543 

8.543 

8.543 

16.51 

24.48 

24.48 

24.48 

24.48 

12.24 

1.469 

2.937 

2.937 

2.937 

2.937 

2.121 

.7320 

.8337 

.9683 

16.31 

22.59 

22.59 

22.59 

22.59 

11.30 

1.469 

2.937 

2.937 

2.937 

2.937 

2.272 

1.471 

1.675 

1.945 

17.48 

22.59 

22.59 

22.39 

22.59 

11.30 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

CRT  MATRIX 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

3.565 

3.565 

3.365 

3.565 

3.565 

199.9 

199.9 

199.9 

199.9 

3.635 

3.565 

3.565 

3.565 

3.565 

.0000 

4.248 

4.248 

4.248 

4.248 

4.248 

12.17 

12.17 

12.17 

12.17 

4.248 

4.248 

4.248 

4.248 

4.248 

.0000 

Figure  C4, 

(Sheet 

3  of 
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VALUE  OF  THE  SRA01ENT  AT  H(J,8I  =  .4583 

VALUE  OF  THE  H(3,8I  *  6.4158 


Figure  C.4.  (Sheet  S  of  5) 
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PROGRAM  LEVEECOR 


C  Corrections  8-12-87 

C  .  . . . . . . . 

C  THE  FOLLONINS  PROGRAM  IS  WRITTEN  BYi 

C  -DR.  THOMAS  F.  WOLFF,  ASSISTANT  PROFESSOR 

C  AMO 

C  -HASDAI  N.  HAJI,  PH.D  CANDIDATE 

C 

c  - 

C 

C  DEVELOPED  AT  THE  A.  H.  CASE  CENTER 

C  COMPUTER  AIDED  EN81NEER1N8  FACILITY  AT 

C  MICHIGAN  STATE  UNIVERSITY  .  SPRIN8  1987 

C 

c  - 

C 

C 

C  THE  PURPOSE  OF  THI8  PROGRAM  IS  TO  ANALYZE  GROUND  NATER  FLOW  BENEATH 
C  FLOOD  CONTROL  LEVEES.  THE  PROGRAM  SOLVES  FOR  THE  DISTRIBUTION  OF  HEAD 
C  FOR  A  THREE  DIMENSIONAL  POROUS  MEDIA  RON  DESCRIBED  BY  LAPLACE  EQUATION. 
C  THE  FINITE  DIFFERENCE  METHOD  IS  USED  -AS  A  MATHEMATICAL  TOOL  FOR  A 
C  A  BY  15  6RI0  MESH  TO  REACH  THE  DESIRED  SOLUTION. 

C  ...  .  . . 


CHARACTER*64  FILEIN 
CHARACTERS  FILEQUT 
CHARACTER»80  TITLE 
1NTE6ER  TRIES, ITT ,R, RN,C,CM 
INTEGER  DONE 
REALM  THETA 
REAL  GRADIENT,!! 

REAL  TEHP,TOPHN,BOTHN 

REAL  L1A,LIB,L2, HR, HL,KF,KBL,KBR, RADIUS, D.COEFF 
REAL  Z, DILEFT tOi 7 ,Oi 16) 

REAL  DIRIGHT(Oi7,Oil6) 

REAL  DYTOP(Oi7,Oi 16) 

REAL  DYBGT(0i7,0i 16) ,DICTR!0i7,0i 16) 

REAL  DYCTR(0i7,0il6) 

DOUBLE  PRECISION  CLTI0i7,0i16),CD(0j7,0i16) 

DOUBLE  PRECISION  CT(0i7,0i16),CRT(0i7,0i16),CB(0i7,0i16I 
DOUBLE  PRECISION  HNfOi7,Oil6),H<Oi7,Oil6),HDIO:7),HO 
REAL  MR,MRN,MRI,HRA,TOL 

C  .  . . . . . 

c 

c  VARIABLE  DICTIONARY 
C  . . 

c 

c 

C  C  INCREMENTAL  INDEX 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


C8  STORES  INTERMID1ATE  VALUE  TO  PE  USED  IN  (IAIN  E6U. 

CO  STORES  1NTERMI01ATE  VALUE  TO  BE  USED  IN  MAIN  EBU. 

CLT  STORES  1NTERMID1ATE  VALUE  TO  BE  USED  IN  (IAIN  EBU. 

CRT  STORES  INTERMEDIATE  VALUE  TO  BE  USED  IN  MIAN  EBU. 

CT  STORES  INTERMEDIATE  VALUE  TO  BE  USED  IN  MIAN  EBU. 

D  ABUIFER  THICKNESS 

DISTANCE  BETWEEN  CONSECUTIVE  COLUMNS  AT  THE  BOTTOM 
OF  THE  MATRIX. 

DTCTR  AVERASE  HORIZENTAL  DIST.  BETWEEN  TWO  CONSECUTIVE  COL. 

DXLEFT  THE  HORIZANTAL  DISTANCE  AT  THE  LEFT  OF  A  BRIO  MESH. 

OtRISHT  THE  HORIZANTAL  DISTANCE  AT  THE  RSHT  OF  A  BRIO  MESH 

DYCTR  AVERABE  VERTICAL  DIST.  BETWEEN  TWO  CONSECUTIVE  RONS. 

DYBOT  VERTICAL  DISTANCE  FROM  ONE  NODE  TO  THE  NETT  LONER  NODE 

DYTOP  VERTICAL  DISTANCE  FROM  ONE  NODE  TO  THE  NEXT  UPPER  NODE 

HR  INITIAL  HEAD  OF  THE  RIVER 

HL  INITIAL  HEAD  BEYOND  THE  LEVEE  (LANDWARD) 

HN  STORES  NEN  HEAD  OF  EACH  NODE  PRODUCED  BY  THE  LAST 
ITERATION 

HO  STORES  OLO  HEAD  VALUES  OF  EACH  NODE  PRODUCES  BY 
PREVIOUS  ITERATION. 

ITT  ITERATION  NUMBER 

KBL  HYORULIC  CONDUCTIVITY 

KBR  HYDRULIC  CONDUCTIVITY 

KF  HYDRULIC  CONDUCTIVITY 

L1A  DISTANCE  BETNEEN  THE  LEFT  SIDE  OF  THE  LEVEE  AND  THE 

RIVER 

LIB  DISTANCE  BETNEEN  THE  RI6HT  SIDE  OF  THE  LEVEE  AND  THE 

RIVER  ! 

L2  NIOTH  OF  THE  LEVEE 

HR  RESIDUAL  HEAD 

MRN  THE  NEN  RESIDUAL  HEAD 

RADIUS  RADIUS  OF  THE  CURVATURE 

THETA  AN6LE  BETNEEN  THE  RIVER  AND  THE  LEVEE  CURVATURE 
TRIES  NUMBER  OF  TRIES 

I  HYDRULIC  QRAOIENT 

Z  THICKNESS  OF  SOIL  LAYERS  I 


C  PROGRAM  HEADIN6 


NRITE  (»,»)  'PROGRAM  LEVEECOR' 

WRITE  (♦,»)  'SEEPAGE  ANALYSIS  AT  LEVEE  CORNERS' 
NRITE  (♦,*) 


C  OPEN  DATA  FILES 


NRITE  (»,»)  'ENTER  INPUT  FILE  NAME' 
READ  (*, ' (A) ')  FIIEIN 
OPEN  I20,FILE=FILEIN) 

REMIND  (201 
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WRITE  (*,*)  'ENTER  OUTPUT  FILE  NAME' 
R£AO<*,'(A)'l  FllEQUT 
OPEN  I 10.FILE-FILEQUT) 

REMIND  (10) 

C  . .  . . 

c 

C  INPUT  PHASE  OF  THE  PROGRAM 
C 


WRITE (10,11) 

WRITE<*,11) 

11  FORMAT !1H1,‘  INPUT  PARAMETERS  OF  THE  PROGRAM* I 
WRITEIlOf 12) 

WRITE!* , 12) 

12  FORMAT (1H0,' - ',/) 

REAS  <20, ‘(A)')  TITLE 

WRITE<*,*1  TITLE 
WRITE  (10,*)  TITLE 
REAO  <20, «)  HR 
MRITE<10,21)  HR 
NR1TE(»,21)  HR 

21  FORMATttHO,'  HR  VALUE* '  ,F9.3,) 

REAO  (20,*)  HL 

WRITEdO ,22)  HL 
WRITE!*, 22)  HL 

22  FORMAT <IH0, 'HL  VALUE*'  ,F9.3,1 
REAO  (20,*)  L1A 

WRITEdO, 23)  LIA 
WRITE!*, 23)  LIA 

23  FORMATdHO,' LIA  VALUES'  ,F9.3,) 

REAO  (20, «)  LIB 

WRITEdO, 24)  LIB 
WRITE!*, 24)  LIB 

24  FORMATI1HO,'  LIB  VALUES'  ,F9.3,) 

REAO  (20, «)  L2 

WRITEdO, 25)  L2 
WRITE!*, 25)  L2 

25  FORMATdHO,'  L2  VALUES'  ,F9.3,) 

REAO  (20,*)  THETA 

WRITEdO, 26)  THETA 
WRITE!*, 26)  THETA 

26  FORMAT !1H0,'  THETA  VALUE) '  ,F9.3,) 

REAO  (20, «)  KF 

WRITEdO, 27)  KF 
WR1TE(*,27) 

27  FORMATdHO,'  KF  VALUES'  ,FV.6,) 

READ  (20,*)  KBL 

WRITEdO, 28)  KBL 
WRITE!*, 28)  KBL 

28  FORMATdHO,'  KBL  VALUES'  ,F9.6,) 

REAO  (20,*)  KBR 
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WRITE ( 10 , 29)  KBR 

NR1TE<*,29)  KBR 

29 

FQRNATUHO,'  KBR  VALUE:"  ,F9.4,) 

READ  (20, O  RADIUS 

NR] TE (10, 301  RADIUS 

MRlTE(t,30)  RADIUS 

30 

FORMAT (IHO,"  RADluS  VALUE:"  ,F9.3,) 

READ  (20,0  D 

NR!TE(10,31)  D 

NRITE(t,3I)  D 

31 

FORMAT (1H0,‘  D  VALUE:'  ,F9.3,) 

READ  (20,0  TOL 

NRITE<10,32)  TOL 

NRIT£(*,32)  TOL 

32 

FORMAT (1H0,'  TOL  VALUE: '  ,F9.4,I 

READ  <20,0  TRIES 

NRITE(I0,33)  TRIES 

NRITE(#,33)  TRIES 

33 

FORMAT (1H0,‘  TOTAL  NUMBER  OF  TRIES:', IS,) 

READ  (20,0  Z 

NRITE  (10,34)  Z 

NRITE<«,34>  Z 

34 

C _ 

FORMAT (1 HO,'  THICKNESS  VALUE:'  ,F12.6,) 

C  INITIALIZATION  PHASE  OF  THE  PROGRAM 
C 


C  CONVERT  THETA  TO  RADIANS 
THETA=THETA/57. 29578 

C  INITIALIZE  HEADS 

DO  40  R=0 , 7 
H(R,0)*0 
H(R,16)=0 

IF  (R.LE.3)  THEN 
HD(R)=HL 
ELSE 

HD(R)«HR 

ENDIF 


CONTINUE 

DO  44,  C=1 , 15 
H  (7 ,0  =0 
H(6,C)*HR 
HO,C)«HL 
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H(0,CI=0 
DO  42,  R=2,5 
H(R,C)=(HR+HL)/2 
42  CONTINUE 

44  CONTINUE 


C  IMTIALIZE  DXLEFT 

W  46,R=1 ,6 
DJLEFT  4Rf 1 )=0 
DO  45  C-2,5 
DXLEFT (R,C) ^RADIUS 

45  CONTINUE 

46  CONTINUE 

X3=SORT ( (KF/KBL) *Z*D> 

IF  (RADIUS.GT.X3)  THEN 
DXLEFT ( t ,&)  «  RADIU3-D«TAN (THETA/2) 

ELSE 

DXLEFT (1,6)  ■  RADIUS- (RADIUS-X3) *TAN (THETA/2) 

ENDIF 

DO  47,R=2,6 
DXLEFT (R,6)=RA0IUS 

47  CONTINUE 

DO  50.RM.2 
DO  49,C=7,10 
DXLEFT  <RfC) =20 
4?  CONTINUE 
50  CONTINUE 

DO  55,C=7 ,10 

DXLEFT (3,CI =RADIUS* (THETA/4) 

55  CONTINUE 

DO  60,C=7,10 

DXLEFT (4,C)=( (RADIUS+L2) * (THETA/4 ) ) 

60  CONTINUE 

DXLEFT (5,7) S(RAD1US*L2*0.5»L1 A*. 125* (L1B-L1A))* (THETA/4) 
OXLEFT(5,8)a(RAOIUS+L2*0.5*LlAL25»<LlB-LI  AIM  (THETA/4) 
DXLEFT(5,9)= (RA0IUS+L2+0. 5*L 1 A+. 375* (L 1B-L 1  A) ) * (THETA/ 4 ) 
DXLEFT (5, 10) =(RADIUS+L2*0,5«L1B)*( THETA/4) 

DXLEFT (6,7)* (RADIUS+L2+L1A+. 25*(L1B-L1 Al )* (THETA/4) 
DXLEFT (6,8)a(RADlUS+L2*LlAL5ML18-LlA)  IMTHETA/4) 

DXLEFT (4,9)  =  (RADIUS*L2*L1A+. 75* (L 1 B-L 1  A) ) « (THETA/ 4 > 
0XLEFK6, 1 01 :  (RAO  (US«L2*L1BI  *(  THETA/4) 
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DO  85 ,R=2, 6 
DXLEFT  IR,11)=RAD1US 
85  CONTINUE 

DO  90,R=1,6 
DO  89tC-12t 15 
DILEFT (R,C)=RADIUS 

89  CONTINUE 

90  CONTINUE 

C  INITIALIZE  DXRISHT 

DO  105,R=1,6 
00  100,C«14tl,-l 

DXRISHT (R,C)aDXLEFT (R,16-C) 

100  CONTINUE 
105  CONTINUE 

DO  U0,R=1,6 
DXRISHT(R,15)*0 
110  CONTINUE 

C _ 

c 

C  THIS  PART  OF  THE  PR06RAN  COMPUTES  DYTQP 

C . . . . 


DO  115,01,5 
0YT0P<1,CI*XJ-RADIUS 
IF  <DYTQP(1,C1.LT.1.0)  THEN 
DYTOP  ( 1 ,0  =D 
ENDIF 
115  CONTINUE 


DO  120,C*ll,15 
DYTOP  ( 1 ,0  =X  J-RAD IUS 
IF  (DYTOPU.O.LT.l.O)  THEN 
DYTQP  (1,0=0 
ENDIF 
120  CONTINUE 

DO  125,0=6,10 

DYTOP  I  l,Cf«DYTOPU,  11  l/COSI  THE!  A/2) 
125  CONTINUE 

DO  126,C=1, 15 

DYTQP(2,CI=RADIUS 


Figure  C5.  (Sheet  6  of  12) 
Cl  7 

I _ 


DYTOP  (3,0=12 

126  CONTINUE 

00  135,  R=4,  5 
DO  127 ,C-1 ,6 
0YT0PlR,C)=LlA/2 

127  CONTINUE 

DYTOP (R, 7J  =L1A  ♦  (0.25ML16-L1AI/2) 

DYTQP(R,B»=L1A  *  <0.5«(L1B-L1A)/21 
DYTQP(R,?)*UA  ♦  (Q.73*(LlB-LlAI/2) 

DO  134,0=10,15 
DYTOP(R,C) =LlB/2 

134  CONTINUE 

135  CONTINUE 

DO  140,0=1,15 
DYT0P(6,d*0 
140  CONTINUE 

0 _ 

C 

0  THIS  PART  OF  THE  PROGRAH  CALCULATES  THE  DYBOT  VALUES 
C  . .  . 


DO  150,0*1,15 
DYBOT  (1,0-0 
150  CONTINUE 

DO  155,R=2,6 
DO  152,C=1,15 
DYBOT  (R,CI=DYTOP(R-1,0 
152  CONTINUE 
155  CONTINUE 

C _ 

c 

c  THE  F0LL0N1NB  PART  CALCULATES  DXCTR  AND  DYCTR 
C 


DO  160,0=1,15 

DICTR(3,0  =  (DXlEFT(3,CHDlRlGHT(3,d)/2 
DYCTR(3,C)=DYB0T(3,O/2 
160  CONTINUE 

DO  170,0=1,15 

DYCTR  (4,0  =  (DYLEFT  (4,C)+DXRIGHT(4,CI)/2 
DYCTR (4, CI=DYTOP  14,0 12 
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170  CONTINUE 


DO  180,C*1,15 
DO  175,11*1,2 

DXCTR(R,C)* (DXLEFT (R ,0 +DXRI BHT ( R , C) ) /2 
DYCTR(R,C)*(DYT0P<R,CHDY8QT(R,C))/2 
175  CONTINUE 
180  CONTINUE 

DO  190,C*i,15 
DO  185,R*5,6 

DXCTR(R,C) * (DXLEFT (R,C) +DXRI6HT (R,C) > /2 
DYCTR (Rf C) * (DYTOP (RtC) +DYBOT (Rf C) ) /2 
185  CONTINUE 
190  CONTINUE 

C 

C 

C  PRINT  MATRICES 
C 


NR1TEU0,191) 

191  FORNAT UH1,' INITIAL  HATRID 

WRITE  110,192) 

192  FORH AT (1  HO,' *******  *****«') 

CALL  DECHO  HAT (H) 

WRITEdO,  193) 

193  FORNAT UHl, 'DXLEFT  MATRIX ' ) 

CALL  ECHO  NAT (DXLEFT) 
NRITE(10,194) 

194  FORMAT ( 1H1 , '0XRI6HT  MATRIX ' > 

CALL  ECHO.NAT(DXRIBHT) 
MR1TE(10,195) 

195  FORMAT  1 1HI,'  DYTOP  MATRIX ' ) 

CALL  ECHO.NAT (DYTOP) 

WR1TE< 10,196) 

196  FORMAT ( 1HI , 'DYBOT  MATRIX ' ) 

CALL  ECHO.NAT(DYBOT) 

WRITE (10,197) 

197  FORNAT (1H1, 'DYCTR  MATRIX') 

CALL  ECHO.NAT (DYCTR) 
NRITE(10,199) 

199  FORMAT ( 1HI , 'OXCTR  MATRIX'! 

CALL  ECHO.NAT (DXCTR) 

C  ... 

c 

C  CALCULATE  C  COEFFICIENTS 

C 

C 
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DO  240,  02,14 
DO  230,02,5 

CLT  4R, 1 ) =0 

CLT (R|C)«D*KF« (DYTOP (R,C>+DYBOT (R,C>  > / (OXLEFT (R,C) *2) 
CLnR,15)*D«KF»<DYTOP<R,15)+DYBOT<R,15))/<0XLEFT(R,l5)*2> 

CT (R , 1 )*D*KF*0XRI6HT (R ,  1 )  / COYTOP (R , 1 ) *21 

CT  <R  ,0  =D*KF*<DXLEFT  <R,C)^0XRI6HT  (R,C) )  /  (DYTOP  (R, 0*2) 

Cl IR , 15) SD*KF*DXLEFT (R , 15) / (DYTOP (R , 15) *2) 

CRT  4R,  1)  »D*KF»  (DYTOP  (R,  D+DYBQT  <R ,  1 ) )  /  (SXRI6HT  (R ,  1  )*2) 

CRT (R,C)*D»KF« <DYTOP(R,C)+DYBQT (R,C) ) / (DXRI8HT (R,C)«2) 
CRT(R,15)*0 

CB (R, 1 ) «D«KF»DXRI8HT (R, 1 ) / (DYBOT (R, 1 >  «2) 
C81R,C)a0*KFHDXLEFT(R,C)+DXRI8HT(R,C))/{DY8QTIR,C)*2) 

CB  I R ,  1 5)  a0*KFOXLEFT  ( R , C)  /  ( DYBOT  1R  ,C )  *2 ) 

IF  (R.LE.3)  THEN 
COEFF*KBL 
ELSE 

C0EFF*KBR 

ENDIF 

CD  1R ,  1 ) a <COEFF«OXCTR  1R,  1 HDYCTR  ( R ,  1 ) )  /  Z 
CD(R,C)a(COEFF*OXCTR(R,C)«DYCTR(R,C))/Z 
CD1R,15)»(C0EFF#DXCTR(R,15)*DYCTR(R,15))/Z 


230  CONTINUE 
240  CONTINUE 

NRITE(10,224) 

224  FORMAT ( iHl , 'CLT  MATRIX') 
CALL  DECHO.NAT(CLT) 

WRITE (10,225) 

225  FORMAT ( IHl , 'CT  MATRIX ' ) 
CALL  DECHO.NAT(CT) 

NRITEU0.226) 

226  FORMAT < IHl t 'CRT  MATRIX ‘ > 
CALL  DECHO.HATICRT) 

WRITE! 10,227) 

227  FORMAT ( IHl, 'CB  MATRIX') 
CALL  DECHO.HATtCB) 

WRITE ( 10 , 228) 
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228  FORMAT (1  HI, 'CO  MATRIX') 

CALL  OECHQ  MAT (CO) 


ITT=0 

D0NE=0 


C  _ 

c" 

C  SET  MEN  HEADS  EQUAL  TO  OLD  HEADS  AT  EDSES  OF  MATRIX 
C  ..  _ 


DO  247,C*l,15 
HN(1,C)»HU,C) 

HH(6,C)*H(6|C) 

247  CONTINUE 

DO  24B,R*lf6 
HN(R,1)*H(R(1) 

HNIR,15)=HIR,15) 

248  CONTINUE 

C 

C  MAIN  ITERATION  LOOP 
C _ 

2481  CONTINUE 

DO  250,R=2,5 
DO  249,C«l,15 

IF  (R.EQ.2.AND.C.8E.7.AND.C.LE.10)  THEN 
HN(R,C)  *  HNI2.A) 

ELSE 

TEMP  *  CLT<R,C)*H<R,C-i)*  CT(R,C)»H(RM,C>  ♦  CRT(R,C)»H(R,C+1) 
TOPHN  »  TEMP  ♦  CD(R,C)*HD(R)  ♦  CB(R,CMH(R-lfC) 

BOTHN  «  CLTIR.C)  ♦  CT(R,C)  ♦  C8(R,C>*  CD(R,C)  ♦  CRTIR.C) 
HN(R,C)  «  TOPHN/BOTHN 
END  IF 

249  CONTINUE 

250  CONTINUE 

c - 

C  RESET  RON  1  HEADS 

C - 

H0-H(3,I)-HL 
DO  251,C=l,15 

HN(l,C)=HL*H0«EXP(-I«(0TT0P(l,U«DyT0P(2tll)/X3) 

251  CONTINUE 
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c _ 

c 

C  CHECK  FOR  HAXlflUn  RESIDUAL 
C 


«R=0.0 

OQ  255,R=l,6 
00  253,C=2,M 
HRN=(HN(R1C)-H(R,CI ) 
MRA=ABS(NRN) 

IF  (MRA.GT.HR)  THEN 
MRX-HRN 
HR-NRA 
RN-R 
CN-C 
ENOIF 
253  CONTINUE 
255  CONTINUE 


C _ 

c 

C  CHECK  FOR  END  DUE  TO  TOLERANCE  OR  TRIES 
C 


IF  ((KR.LT.TOl).OR. (1TT.6T.TRIES))  THEN 
NRITEU0.270) 

270  FORITATUHl, ‘FINAL  H  MATRIX') 

CALL  DECHO.HATIH) 

MR1TE(10,275)  ITT 
NRITE(»,275I  ITT 

275  FORMAT (IH0,‘  TOTAL  NUMBER  OF  ITERATIONS  ',110) 

NRITEOO,*)  TITLE 
MR1TE0,*)  TITLE 
NRITEI*,*) 

6RA0IENT«H(3|8)/Z 
WRITE ( 10 ,27&>  GRADIENT 
MRITE  0,276)  GRADIENT 

276  FORMAT  1 1 HO VALUE  OF  THE  GRADIENT  AT  H(3,8)  *  ',F9.4) 

MRITE(10,277)  H(3,8) 

NRITEO,277)  H(3,8) 

277  FORMAT ( 1  HO, ‘  VALUE  OF  THE  H(3,8)  «  ' ,F9.4) 

REMIND  (20) 

CLOSE  (20) 

REMIND  (10) 

CLOSE  (10) 
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B0NE=1 


ELSE 

00  400 ,R=i ,5 
00  300,01,15 
H(R,C)=HN(R,C) 

300  CONTINUE 

400  CONTINUE 

ITTMTTd 

NRITE(*,4101  ITT,  H(3,8),  HRX ,  RH,  CH 
410  FORMAT  (IX, '  ITERATION  ',15, 

V  '  HEAD  «  * ,F9. 5, *  MAX  RES  *  ',F9.5,'  AT  ',13, ',‘,131 


ENOIF 

IF  (DONE  .EQ.  0)  GOTO  2481 

STOP 

END 


C 

C  THIS  SUBROUTINE  DISPLAYS  THE  CONTENT  OF  A  MATRIX. 
C 

SUBROUTINE  ECHO  HAT (HAT) 

REAL  MAT (Oi 7 ,0i 16) 

INTEGER  I,J 

DO  10,1=1,6 

WRITE (10,15)  (HAT(I,JI, 01,15) 

15  FORHATdX,  15(610.41) 

10  CONTINUE 
RETURN 
END 

SUBROUTINE  DECHO.HATIHAT) 

DOUBLE  PRECISION  NAT(0:7,0il&) 

INTEGER  1,J 

DO  10,1=1,6 

WRITE! 10, IS)  (HAT ( I , J ) , J=l , 15) 

15  FORHATdX,  15(610.4) ) 

10  CONTINUE 
RETURN 
END 
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APPENDIX  D:  COMPARISON  OF  PROGRAMS 


Example  Problem 


1.  In  this  appendix,  a  comparison  is  made  among  the  solutions  obtained 
for  a  common  example  problem  using  the  computer  programs  LEVEE3L,  LEVEEIRR, 
LEVEECOR,  and  the  conventional  method.  The  example  selected  is  the  same  as 
that  used  for  Figure  45  in  the  main  text.  Foundation  conditions  are  assumed 
to  consist  of  a  10-ft-thick  top  blanket  overlying  an  80-ft-thick  substratum. 
The  ratio  of  the  horizontal  substratum  permeability  to  the  vertical  top 
blanket  permeability  is  assumed  to  be  1,000.  The  foundation  is  assumed  to 
extend  1,500  ft  riverside  of  the  levee  to  an  open  seepage  entrance.  The  levee 
base  width  is  assumed  to  be  300  ft.  The  foundation  is  assumed  to  extend 
infinitely  in  the  landside  direction.  The  calculated  residual  head  at  the 
levee  toe  using  the  conventional  method  is  8.97  ft. 


Comparative  Solutions 


2.  The  heads  obtained  using  the  computer  analyses  and  conventional 
analyses  are  plotted  as  a  function  of  horizontal  distance  in  Figure  Dl. 
Computer  solutions  for  the  residual  head  at  the  levee  toe  range  from  6.6  to 
8.9  ft,  a  difference  of  2.3  ft.  Two  different  segment  assumptions  are  shown 
for  LEVEEIRR  in  Figure  Dl.  Computer  solutions  in  this  example  reflect 
generally  lower  heads  than  the  conventional  solution,  a  trend  also  noted  in 
the  case  studies  cited  in  the  main  report.  These  variations  are  larger  than 
was  generally  expected  at  the  outset  of  the  research.  Reasons  for  these 
variations  are  not  completely  understood  at  the  time  of  preparing  this  report; 
there  are,  a  number  of  factors  that  may  contribute.  These  are  discussed 
below. 


Differences  in  Program  Assumptions 

3.  The  programs  and  the  conventional  method  employ  certain  fundamen¬ 
tally  different  assumptions  to  derive  the  flow  equations.  In  the  conventional 
method,  flow  is  assumed  to  be  horizontal  in  one  direction  in  the  substratum 
and  is  assumed  to  be  vertical  in  the  top  stratum.  In  LEVEE3L,  both 
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horizontal  and  vertical  flow  is  accounted  for  in  each  of  three  layers.  In 
LEVEEIRR,  flow  assumptions  are  generally  the  same  as  the  conventional  solu¬ 
tion,  but  the  flow  equation  is  still  satisfied  only  approximately  at  a  finite 
number  of  points  rather  than  exactly  at  all  points.  In  LEVEECOR,  flow  in  the 
substratum  is  horizontal  only,  but  has  two  horizontal  components,  and  flow  in 
the  top  stratum  is  vertical  only.  In  LEVEE3L  and  LEVEEIRR,  boundaries  at  an 
infinite  distance  must  be  modeled  by  using  rather  large  finite  distances. 
Results  will  vary  depending  on  the  boundary  distance  specified.  In  LEVEECOR, 
which  most  closely  matched  the  conventional  solution  at  the  levee  toe,  the 
head  on  the  landside  boundary  is  forced  to  be  consistent  with  the  conventional 
closed-form  solution. 


Numerical  Problems  Common  to  All  Programs 


4.  In  each  of  the  three  programs,  the  solution  is  obtained  by  itera¬ 
tion.  As  the  tolerance  is  made  smaller  and  the  number  of  iterations  is 
increased,  a  more  accurate  solution  is  expected.  However,  this  may  not  always 
be  the  case  as  round-off  errors  may  accumulate  where  the  tolerance  or  maximum 
allowable  residual  is  set  extremely  low.  Recommended  tolerance  values  based 
on  limited  experience  with  the  programs  are  given  in  Appendices  A  through  C. 
Associated  with  any  finite  difference  solution  is  a  discretization  or  trunca¬ 
tion  error,  due  to  the  fact  that  a  continuous  system  is  replaced  with  a  finite 
system.  For  the  three  computer  programs  developed,  the  gradients  at  a  node 
point  are  represented  by  a  first-order  approximation  of  the  derivative  of  head 
with  respect  to  distance.  The  truncation  errors  could  conceivably  be  reduced 
and  more  accurate  solutions  obtained  if  higher  order  approximations  were  used 
or  if  node  spacing  was  significantly  reduced.  Both  of  these  approaches  would 
increase  program  complexity  or  execution  time.  The  present  research  was 
viewed  as  an  initial  evaluation  of  the  utility  of  applying  numerical  methods 
to  underseepage  analysis  and  involved  a  deliberate  effort  to  produce  working 
versions  of  three  different  programs  that  were  fundamentally  as  simple  as 
possible . 
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APPENDIX  E:  NOTATION 


c 

d 


H 

*bl 

kbr 

kf 

klh 

kIv 

k2h 

k2v 

k3h 

k3v 

L1 


J3 

q 


z 

z 


1 

2 


A  constant;  inverse  of  the  effective  exit  distance 
Thickness  of  pervious  substratum 
Residual  head  at  levee  toe 

Net  head  at  distance  x  landward  of  levee  toe 
Net  head  on  levee 
Permeability  of  top  blanket 
Permeability  of  landside  top  blanket 
Permeability  of  riverside  top  blanket 
Permeability  of  pervious  substratum 

Horizontal  permeability  of  top  blanket  in  three  layer  system 

Vertical  permeability  of  top  blanket  in  three  layer  system 

Horizontal  permeability  of  middle  stratum  in  three  layer  system 

Vertical  permeability  of  middle  stratum  in  three  layer  system 

Horizontal  permeability  of  substratum  in  three  layer  system 

Vertical  permeability  of  substratum  in  three  layer  system 

Distance  from  riverside  toe  of  levee  to  open  or  blocked  seepage 
entrance 

Distance  from  riverside  toe  of  levee  to  landside  toe 

Distance  from  landside  toe  of  levee  to  open  or  blocked  seepage  exit 

Flow 

Distance  from  landside  levee  toe  to  effective  seepage  entrance 
Distance  from  riverside  toe  of  levee  to  effective  seepage  entrance 

Distance  from  riverside  toe  of  levee  to  landside  toe 

Distance  from  landside  toe  of  levee  to  effective  seepage  exit 
Thickness  of  top  blanket 

Thickness  of  top  blanket  in  three-layer  system 

Thickness  of  middle  stratum 
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